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ABSTRACT 


The  research  goals  were  (1)  to  develop  a  dynamic  response  model  or 
constitutive  relations  for  porous  materials,  (2)  to  test  this  model 
against  results  from  shock  wav<>  experiments,  and  (3)  to  Incorporate  the 
model  into  a  computer  program  lor  the  analysis  of  shock  wave  propagation 
arising  either  from  impacts  or  from  radiation  deposition.  The  model 
developed  exhibits  melting,  vaporization,  and  temperature-dependent 
compaction  resistance,  yielding,  consol  idatio.),  and  spalling.  The  model 
accounts  for  loading  and  unloading,  heating  and  cooling,  and  any  combina¬ 
tion  of  these  processes.  The  computer  program  (SRI  PUFF  1)  for  analyzing 
shock  wave  propagation  problems  employs  artificial  viscosity  in  a  modified 
Lax-Wendroff  integration  scheme. 

Flyer  plate  impact  experiments  (shock  attenuation  tests)  were 
conducted  on  samples  of  porous  iron,  copper,  and  tungsten,  and  stress 
histories  were  recorded  at  the  rear  face  of  the  targets.  Quasi-static 
one-dimensional  compression  tests  were  made  on  samples  of  the  same  porous 
metals,  using  pressures  up  to  10  kbar.  Numerical  values  of  parameters 
in  the  model  were  obtained  from  the  quasi-static  tests  and  from  previously 
reported  Hugoniot  experiments.  Stress  histories  computed  with  the  code 
were  compared  with  stress  records  obtained  from  the  attenuation  experi¬ 
ments.  The  computed  (predicted)  peak  stresses  and  arrival  times  of  the 
waves  generally  agreed  to  within  20  percent  or  better  with  the  measured 
values.  Precursor  amplitudes  and  other  wave  front  features  were  correctly 
represented  in  the  computed  histories.  The  model  appears  to  adequately 
represent  the  dynamic  response  of  porous  materials,  the  major  uncertainty 
being  the  numerical  values  of  parameters  for  a  given  material. 
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SECTION  I 
INTRODUCTION 


Previous  work  (Refs.  1-5)  performed  for  the  Air  Force  Weapons 
Laboratory  (AFWL)  has  demonstrated  that  porous  solids  (also  called  dis¬ 
tended  solids,  or  foams)  can  be  very  effective  in  attenuating  the  peak 
stress  of  a  propagating  stress  pulse.  Hence,  to  protect  a  structural 
material  from  shock  damage  caused  by  short-duration  impulsive  loadings, 
an  external  layer  of  a  porous  material  might  be  used.  For  effective  use 
of  such  porous  materials,  it  is  necessary  to  be  able  to  predict  the 
response  of  the  material  to  expected  shock  loadings.  This  report  provides 
a  procedure  for  predicting  the  response. 

Research  performed  at  Stanford  Research  Institute  (SRI)  and  else¬ 
where  (see  Refs.  1-20,  for  example)  has  provided  a  considerable  body  of 
Hugoniot  and  other  dynamic  data  on  porous  metals  and,  to  some  extent, 
on  porous  plastics  and  ceramics.  To  make  effective  use  of  such  data  in 
vulnerability  studies  and  survivability  design  of  possible  reentry 
systems,  it  is  necessary  to  employ  this  data  to  construct  a  mathematical 
model  of  the  dynamic  response  of  porous  materials  to  impulsive  loading. 

To  make  the  needed  predictions,  the  model  must  be  incorporated  into  a 
computer  code  (such  as  PUFF  66)  for  analyzing  shock-wave  propagation. 

The  present  effort  was  initiated  to  develop  a  prediction  capability  for 
shock  propagation  and  attenuation  behavior  in  porous  solids  using  data 
(from  AFWL  TR-68-335)  on  porous  copper,  iron,  and  tungsten.  This  goal 
has  been  met  with  development  of  the  SRI  PUFF  1  computer  code. 


*  As  used  in  this  report  the  term  "Hugoniot  is  applied  to  foams  to 
mean  the  locus  of  final  macroscopic  pressure-volume  or  pressure- 
particle  velocity  states  of  shocked  material  deduced  from  experimental 
observations  and  invocation  of  the  so-called  "jump"'  conditions  'vhich 
express  conservation  of  mass  and  momentum  across  a  shock  fror*.  This 
Hugoniot  represents  actual  rather  than  average  states  of  the  material 
only  if  the  scale  of  porosity  is  such  that  equilibrium  is  attained 
within  the  time  scale  of  the  experiment. 


In  addition  to  experiments  on  porous  metals,  three  exploratory 
dynamic  experiments  were  performed  at  peak  stresses  below  50  kbar  on 
"3-D"  quartz  phenolic  material  supplied  by  AFWL.  Additional  experiments 
are  currently  being  performed.  To  present  the  results  of  these  experi¬ 
ments  in  the  most  useful  manner,  the  results  of  the  three  quartz 
phenolic  experiments  performed  under  the  present  contract  will  be 
presented  in  a  later  report. 


SECTION  II 
BACKGROUND 

1 .  THEORETICAL 

The  historical  development  of  equations  of  state  for  porous 
materials  is  closely  connected  with  the  development  of  computer  programs 
for  wave-propagation  calculations.  In  this  section,  we  outline  the  develop 
ment  of  wave-propagation  computer  programs  and  then  describe  in  detail 
equation-of-state  formulations  that  have  been  developed. 

In  1950,  von  Neumann  and  Richtmyer  (Ref.  21)  initiated  the  artificial 
viscosity  (or  Q)  method  for  using  digital  computer  codes  to  solve  the 
equations  of  wave  propagation.  With  this  technique,  infinitely  steep 
shock  fronts  are  not  allowed  to  develop  and  the  entire  field  can  be 

treated  as  a  continuous  flow.  Shock  fronts  appear  as  regions  of  high 
stress  gradient  but  not  as  discontinuities.  The  artificial  viscosity 
tends  to  dampen  all  oscillations  or  perturbations  in  the  flow  field. 

Several  integration  schemes  based  on  the  Q  method  have  been  developed, 
notably  the  Lax-Wendroff  method  (Ref.  22),  the  Runge-Kutta-Gill  method 
(Ref.  ?3),  and  the  "leapfrog''  scheme  (Ref.  21),  which  is  used  in  most 
PUFF  codes. 

The  present  li».e  of  PUFF-type  codes  originated  around  1958  with  the 
development  of  the  SHARK  (Ref.  24)  and  SHARP  (Ref.  25)  codes.  With  later 
developments  at  the  Air  Force  Special  Weapons  Center,  Kirtland  Air  Force 
Base,  the  generic  name  PUFF  was  given  to  the  programs.  Recent  versions 
include  PUFF  (Ref.  26,  27,  28),  PUFF  III  (Ref.  29),  PUFF  IV  (Ref.  30), 

PUFF  IV-EP*  (Ref.  31),  PUFF  V-EP  (Ref.  32),  PUFF  VTS**  (Ref.  33),  FOAM 
PUFF  (Ref.  34),  PUFF  66  (Ref.  35),  and  P  PUFF  66  (Ref.  35).  The  present 
code  (SRI  PUFF  1)  began  as  a  modification  of  the  last  two.  Most  of  the 
PUFF  codes  have  been  described  in  classified  reports,  so  their  charac¬ 
teristics  cannot  be  outlined  here.  A  useful  review  of  the  capabilities 
of  each  of  these  codes  has  been  provided  by  Bothell  and  Archuleta  (Ref.  31) 


*  EP  stands  for  elastic-plastic. 

**  Variable  time  step. 
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Other  PUFF-type  codes  are  available  under  the  names  of  WONDY  (Ref.  36), 

AFOAM  (Ref.  4),  and  AFTON  (Ref.  37).  All  of  the  PUFF-type  codes  use 
artificial  viscosity  and  the  leapfrog  integration  scheme. 

The  attempt  to  insert  into  a  PUFF  code  a  so-called  equation  of  state 

* 

for  a  porous  material  has  in  the  past  met  with  considerable  difficulty. 
First,  the  formulation  of  the  equation  of  state  has  been  hampered  by  a 
lack  of  specific  Information  on  material  behavior.  Most  attempts  to 
develop  so-called  equations  of  state  for  porous  materials  have  been 
directed  toward  plate-impact  problems.  Hence,  radiation  deposition  was 
not  included,  and  it  was  assumed  sufficient  to  provide  only  Hugoniot 
paths  rather  than  the  more  complete  P-V-E  (pressure-volume-energy) 
equation  of  state.  Second,  the  computations  have  been  fraught  with 
oscillations  and  instabilities  associated  with  the  large  compressibilities 
of  the  porous  materials.  This  problem  with  the  computations  points  out 
an  advantage  of  developing  the  computer  program  and  the  equation  of  state 
concurrently. 

The  equations  of  state  postulated  for  porous  materials  have  grown 
progressively  more  complex  as  more  has  been  learned  about  the  materials 
themselves.  At  present,  the  material  behavior  during  shock  compression 
in  plate  impact  experiments  is  the  aspect  best  understood.  The  thermo¬ 
dynamic  stress-volume  paths  followed  by  a  porous  material  during  the 
impact  are  probably  like  that  shown  in  Fig.  la.  In  most  porous  materials 
there  is  an  initial,  approximately  elastic  behavior  up  to  a  somewhat 
ill-defined  "yield  point."  This  yield  point  defines  the  amplitude  of  the 

n  H *  ** 

elastic  precursor  in  the  material.  For  higher  stresses, the  particles 
yield  or  crush, and  pores  are  filled  in  so  that  the  behavior  is  largely 

*  The  pressure-volume-energy  relations  that  describe  the  behavior  of  porous 
materials  do  not  actually  constitute  a  true  equation  of  state  (this 
point  is  discussed  in  detail  later).  However,  because  these  relations 
take  the  place  of  an  equation  of  state,  they  are  popularly  called  an 
equation  of  state. 

**  tfhen  a  load  is  applied  to  a  porous  material,  high  stress  concentrations 
occur  at  contact  points  of  the  particles.  Thus,  even  the  "elastic" 
precursor  will  generally  cause  a  small  amount  of  localized  plastic  flow 
that  enlarges  contact  surfaces  and  relieves  the  high  stress  concentra¬ 
tions  . 
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SOLID  COMPRESSION  PATH 


.COMPRESSION  PATH 
IN  POROUS  MATERIAL 


,UNLOAD  AND 
RELOAD  PATHS 
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OR  YIELD  POINT/ 


SPECIFIC  VOLUME 

(a)  COMPRESSION  PATH  OF  POROUS  MATERIAL 


VOLUME 

(C)  ELASTIC/ 

PERFECTLY  PLASTIC 
RIGID  LOCKING 
MODEL 


VOLUME 

(b)  SIMPLE  LOCKING  OR 
SNOWPLOW  MODEL 


VOLUME 

(d)  ELASTIC/ 

PERFECTLY  PLASTIC 
LOCKING  MODEL 


VOLUME 

(e)  ELASTIC -LOCKING 
MODEL 


FIG.  1  COMPRESSION  PATH  OF  A  POROUS  MATERIAL  AND  FOUR  SIMPLE 
MODELS  FOR  POROUS  MATERIALS 
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irreversible.  After  consolidation  (full  compaction  to  the  solid  state), 
the  compression  path  for  subsequent  compressive  loadings  is  taken  to 
coincide  with  that  of  the  solid  material  at  the  same  temperature.  Simple 
'  that  represent  this  compressive  behavior  in  a  porous  material 
jhown  in  Figs,  lb  through  le. 

Various  models  can  be  used  to  define  Hugoniots,  isentropes,  and  other 
stress-strain  paths,  depending  on  the  experimental  conditions.  The  simple 
locking  or  snowplow  model  specifies  compaction  at  zero  stress  up  to 
consolidation,  and  may  therefore  be  used  if  the  elastic  compression  of 
the  solid  is  unimportant  and  material  compacts  at  an  essentially  conatant 
stress,  or  if  the  shock  stress  is  well  above  th^  consolidation  stress  of 
the  material.  The  elastic/perfectly  plastic  1  :king  i.odel  takes  into 
account  the  elastic  behavior  of  the  solid,  but  does  not  allow  for  any 
strength  of  partially  compacted  material.  The  elastic-locking  model 
does  allow  for  such  strength  and  is  applicable  at  stresses  below  consolida¬ 
tion.  All  of  these  models  lead  to  fairly  straightforward  wave-propagation 
calculations;  therefore,  they  have  been  used  extensively  for  computations 
of  waves  in  porous  materials.  In  the  following  discussion, we  mention  some 
of  the  more  complex  formulations  for  porous  materials. 

In  1965,  Thouvenin  (Ref.  38)  presented  his  plate-gap  model  for  shock 
compression  of  porous  materials  at  high  pressures.  This  model  represents 
a  porous  material  as  a  series  of  plates  separated  by  gaps.  Wave  propaga¬ 
tion  within  the  plates,  impact  of  the  plates  on  each  other,  and  consequent 
removal  of  gaps  were  considered  in  detail  to  determine  the  macroscopic 
behavior  of  the  model.  Thouvenin  considered  particularly  the  impact  of 
a  solid  flyer  plate  on  a  porous  target  made  of  plates  of  the  same  material. 
He  determined  from  wave-propagation  calculations,  that  for  the  same 
impact  velocity,  the  free-surface  velocity  of  a  porous  material  should 
be  equal  to  that  of  a  full-density  solid  of  the  same  material.  He  also 
indicated  that  a  porous  solid  under  the  influence  of  shock  reaches  an 
equilibrium  state  that  is  on  the  same  pressure-volume  Hugoniot  as  that  of 
the  compact  solid,  that  is,  that  energy,  pressure, and  volume  points 
coincide  after  consolidation.  But  the  energy  attained  in  the  material  . 
while  passing  from  the  initial  specific  volume  VD  at  atmospheric  pressure 
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to  the  final  volume  V  is 

E  =  \  S(VQ  -  V)  +  Eq 

where  S  =  stress  of  the  shock  wave 

EQ  =  initial  specific  energy  of  the  material,  usually  negligible 

The  final  energies  behind  a  shock  of  a  given  stress  S  must  therefore 
differ  if  the  initial  volumes  differ.  Hence,  the  coincidence  of  poj ous 
and  solid  pressure-volume  Hugoniot  is  a  violation  of  the  first  law  of 
thermodynamics  and  also  contrasts  with  the  findings  of  experimentalists 
for  strong  shocks.  In  his  analysis, Thouvenin  does  not  take  into 
account  the  difference  between  the  thermodynamic  stresses  and  the 
mechanical  stresses  in  the  wavefront.  If  this  difference  were  considered 
and  the  mechanical  stress  would  follow  the  Rayleigh  line  in  the  P-V 
plane,  then  a  different  Hugoniot  would  be  formed  corresponding  to 
each  initial  density  of  solid  or  porous  material. 

J.  F.  Heyda  (Ref.  39)  modified  Thouvenin's  model  slightly  to 
guarantee  momentum  conservation  across  the  shock  front.  In  this  modified 
model  the  pres sure- volume  Hugoniot  depends  on  initial  density,  unlike 
Thouvenin's  model.  Heyda's  results  correlated  well  with  some  experi¬ 
mental  data  at  high  pressures.  However,  this  model  implicitly  assumes 
an  isentropic  Hugoniot  for  the  solid  plates.  This  may  not  be  serious 
at  low  stresses,  but  becomes  more  in  error  as  the  shock  stress  is 
increased.  Furthermore,  in  a  real  material, with  porosity,  one  normally 
expects  considerably  more  entropy  production  (because  of  microjetting 
and  other  effects)  than  that  provided  by  a  simple  plate-gap  model  even 
when  corrected  for  noncoincidence  of  solid  Hugoniots  and  release  adiabats 

Prindle  (Ref.  40)  also  presented  a  theoretical  study  of  shocks  in 
porous  materials  based  on  the  approach  of  Thouvenin.  He  developed  a 
considerably  more  complex  set  of  equations  to  explain  the  response  of 
Thouvenin's  model  to  compaction  waves.  He  included  material  strength  so 
that  a  precursor  wave  could  be  produced.  Low-pressure  shocks  were 
assumed  to  remove  voids  in  the  direction  of  propagation  but  to  leave 
lateral  voids.  Stronger  shocks  were  assumed  to  squeeze  the  material 
and  remove  the  voids.  The  precursor  wave  was  associated  with  the  first 


type  of  void  removal,  while  the  main  shock  produced  complete  compaction. 
Prlndle  obtained  analytic  solutions  for  Impact  problems.  However,  his 
model  is  based  on  the  same  isentropic  behavior  as  Th-iuvenin's  model,  and 
therefore  suffers  from  the  same  theoretical  inconsistencies. 

Butcher  Ref.  41)  discussed  the  crush-up  or  consolidation  behavior 
for  several  models  that  nave  been  postulated.  He  discussed  in  some 
detail  the  model  that  Fritz  and  Tstylor  (Ref.  42)  developed  for  adiprene 
rubber  foam.  The  model  is  based  on  the  concept  of  a  limiting  density  and 
on  the  static  compression  curve  for  the  material.  The  equation  of  state 
resulting  from  the  analysis  resembles  that  for  a  soft  solid:  it  does  not 
show  the  large  permanent  deformations  that  occur  with  distended  metals 
or  plastics.  Butcher  also  discussed  the  snowplow  model  and  Thouvenin's 
plate-gap  model  and  made  a  few  calculations  with  each  to  illustrate  the 
results  obtainable. 

Models  based  directly  on  Hugoniot  measurements  have  been  proposed 
by  Wagner,  Brooks,  and  Bjork  (Ref.  43),  by  Herrmann  (Refs.  44,  45);  and 
by  Linde  and  Schmidt  (Refs.  3,4),  who  have  discussed  phenomenological 
models  that  have  various  degrees  of  complexity. 

None  of  the  models  discussed  above  explicitly  allow  for  external 
heating,  and  thus  they  are  not  suitable  for  radiation  deposition  problems, 
and  therefore  fail  to  account  completely  for  the  effects  of  internal 
energy.  These  models  may  be  used  for  plate  impact  calculations  in  which 
the  Hugoniot  of  the  solid  matrix  material  and  the  Hugoniot  of  the  consoli¬ 
dated  porous  material  essentially  coincide,  i.e.,  where  thermal  effects 
are  negligible.  Radiation  deposition  computations  cannot  be  made  with 
these  models.  An  early  model  which  included  external  heating  was 
presented  by  Allen  et  al.  (Ref.  34)  in  1964  for  a  material  such  as 
styrofoam.  The  energy-pressure-volume  model  was  incorporated  into 
FOAM  PUFF,  a  computer  program  of  the  PUFF  VTS  type.  The  yield,  or  upper, 
surface  of  the  equation  of  state  for  porous  material  was  similar  to  the 
one  presented  in  this  report.  The  variation  in  the  P-V  plane  was, 
however,  given  by  a  single  exponential  equation.  In  the  direction  of 
increasing  internal  energy,  the  pressure  (or  strength  of  the  foam) 
decreased  linearly  with  energy,  reaching  zero  when  the  internal  energy 


reached  the  sublimation  energy.  Hence,  melting  was  not  specifically 
provided  for.  The  material  was  allowed  to  unload  and  reload  in  h  physi¬ 
cally  reasonable  manner  and  with  a  modulus  like  that  of  the  virgin 
material.  More  recent  models  for  porous  materials  have  been  formulated 
by  other  investigators  for  very  high-pressure  work,  to  emphasize  the 
correspondence  between  equations  of  state  for  porous  and  solid  materials, 
or  to  simplify  the  computations.  However,  the  recent  models  do  not 
seem  to  have  made  important  improvements  to  the  FOAM  PUFF  model  and 
they  are  less  complete. 

Herrmann  (Ref.  44  )  presented  a  formulation  for  an  E-P-V  equation  of 
state  for  crushable  distended  materials.  The  equation  of  state  is  based 
essentially  upon  the  Mie-Grune  sen  equation  of  state  for  solid  materiel 
with  the  addition  of  a  parameter  alpha,  the  ratio  of  the  solid  to  the 
porous  densities.  Herrmann  also  gave  an  analytical  form  for  the  variation 
of  alpha  with  compaction,  thus  providing  a  Hugoniot  for  any  porous  material. 
Herrmann's  complete  formulation  includes  analytical  formulas  for  the 
wave  velocity  during  loading  of  the  virgin  material  and  during  unloading 
and  subsequent  reloading'.  His  sample  calculations  showed  considerable 
erosion  of  the  wave  fronts.  He  noted  that  for  computations  with  porous 
materials,  considerably  longer  running  times  than  for  solids  are  required, 
and  higher  values  of  the  artificial  viscosity  are  needed.  Herrmann 
(Ref.  45)  later  altered  the  formulation  of  the  equation  of  state  to 
provide  for  arbitrary  initial  wave  velocities  of  the  virgin  material. 

In  his  papers .Herrmann  has  specifically  formulated  the  equation  of 
state  only  in  the  P-V  plane:  he  notes  that  the  variation  in  the  energy 
direction  has  not  yet  been  defined.  However,  because  the  porous  equation 
of  state  is  closely  connected  with  the  solid  equation  of  state,  the 
variation  in  the  energy  direction  is  defined  by  default.  In  any  plane 
of  constant  E  in  the  E-P-V  space,  the  P-V  curve  has  the  same  shape  as 
at  E  =  0  (see  Fig.  2).  If  the  porous  material  were  heated  at  constant 
volume  it  would  increase  in  strength,  never  melting,  but  becoming  solid 
at  some  point.  Such  behavior  is  not  physically  reasonable.  Therefore, 
the  Herrmann  equation  of  state  should  be  used  only  for  impacts  where 
internal  energies  less  than  melting  are  of  interest,  or  in  computations 
where  the  stresses  are  well  above  the  consolidation  pressure. 
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FIG.  2  ENERGY-PRESSURE-VOI  UME  (E-P-V)  DEPICTION  OF  HERRMANN’S  EQUATION 
OF  STATE 


As  an  extension  of  Herrmann's  model,  Johnson  (Ref.  46)  proposed  a 
theory  for  strain  rate-dependent  behavior  of  porous  solids.  He  suggested 
that  the  experimentally  observed  slow  rise  of  the  second  wave  in  porous 
materials  results  from  a  significant  amount  of  strain  rate  dependence. 
Johnson's  formulation  starts  with  the  analytical  equation  of  the  Fugoniot 
given  by  Herrmann.  He  then  presents  two  possible  stress  relaxation 
functions  (F  and  F  )  that  are  proportional  to  the  plastic  strain  rate. 

i.  <2 

These  two  functions  are  as  follows: 
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Either  of  these  two  plastic  strain  functions  can  be  used  in  the  equation- 
of-state  relations  as  follows: 


d£  M  dV 
dt  "  V  dt 


(3) 


where 

t  =  time 

F  =  a  strain  rate  function 

M  =  the  modulus  for  the  usual  rate-independent  equation  of  state 
da  *  MdV/V. 

Johnson  calculates  stready-state  wave  profiles  in  materials  for  these 
two  stress  relaxation  functions.  For  both  functions  the  thickness  of 
the  shock  front  decreases  markedly  with  increased  stress.  The  profiles 
are  similar  for  the  two  formulations  but  the  one  that  is  a  function  of 
stress  seems  to  be  steeper  at  high  stresses. 

Of  the  many  models  for  quasistatic  behavior  of  porous  solids,  one 
empirical  equation  is  worthy  of  mention  here.  Kawakita  and  Tsutsumi 
(Ref.  47)  presented  an  empirical  equation  of  state  for  powder  compression. 
The  equation  is  for  static  work  and  corresponds  to  an  isotherm  rather 
than  to  the  Hugoniot: 


K 


V  -V 
o 


a  b  P 
1  b  P 


(4) 
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where  K  Is  the  bulk  modulus  and  a  and  b  are  constants.  The  equation 
was  fitted  to  data  on  metal  powders,  ferrite  powders,  foodstuffs,  and 
other  powders,  and  appeared,  in  general,  to  give  a  good  representation 
of  the  data. 

Biot  (Ref.  48)  developed  equations  of  state  for  elastic  and  visco¬ 
elastic  behavior  of  fluid-filled  anisotropic  porous  media.  However, 
the  elastic  moduli  of  the  composite  material  were  not  given  as  functions 
of  the  known  parameters. 


Mackenzie  (Ref.  49)  derived  expressions  for  the  elastic  constants 
of  a  solid  containing  spherical  holes.  It  was  intended  that  the  results 
would  be  applicable  to  sintered  material  because  on  sintering  the  pores 
usually  assume  an  approximately  spherical  shape.  The  analysis  was  made 
for  static  elastic  behavior  only.  The  results  of  the  calculations  were 
formulas  for  bulk  and  shear  moduli,  which  increase  with  density  of  the 
material,  reaching  values  for  the  solid  when  the  density  is  that  of  a 
solid.  The  equations  for  the  moduli  are  as  follows: 


1  P8  3(1  -  PSPs'  P, 

K  "  K  p  +  4  G  P 

0  0 


(5a) 


G  -  G 

-  =  5(1  -  P/p  ) 

e  s 

o 

where 

K  and  G  *  the  bulk  and  shear  moduli  of  the  porous  material 

K  and  G  =  the  moduli  of  the  solid  material 
o  o 

p  and  p  =  densities  of  the  porous  and  solid  materials,  respectively, 
s 

These  formulas  show  a  linear  variation  of  the  shear  modulus  with  density. 
The  bulk  modulus  increases  slowly  initially  during  compaction  and  then 
rapidly  as  the  density  approaches  that  of  the  solid. 
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2.  EXPERIMENTAL 

A  considerable  body  of  experimental  shock  propagation  data  for 
porous  materials  has  been  obtuined.  Selected  important  experimental 

investigations  are  summarized  below. 

a .  Soil  and  Rock 

Grogson,  Ahrens,  and  Petersen  (Ref.  6)  studied  quartzite,  sandstone, 
basalt,  calcite,  and  plagioclase  rocks  in  the  range  from  4  to  250  kbar. 

High  values  of  the  Hugoniot  elastic  limit  were  observed  in  solid  rocks  - 
50  to  100  kbar  in  quartzite,  40  to  50  kbar  in  feldspar  and  basalt,  and  20 
kbar  in  calcite  and  marble.  Limestone  (a  porous  rock  initially  at  87  to 
95  percent  of  solid  density)  showed  elastic  precursors  of  10  to  15  kbar  and 
compaction  to  the  solid  (calcite)  pressure-volume  Hugoniot  at  25  to  40  kbar. 
The  sandstones  (74  to  88  percent  of  solid  density)  tested  exhibited  pre¬ 
cursors  of  10  kbar  and  compaction  to  the  solid  (quartz)  Hugoniot  at  50  to 
100  kbar. 

/nderson  et  al.  (Ref.  2)  conducted  a  study  of  equations  of  state  of 
poroi  i  earth  media.  The  study  included  experiments  on  wet  playa,  with 
ini*  ;.sl  dry  densities  of  58  and  73  percent  of  solid  density,  a  discussion 
of  several  forms  of  the  Mie-Gruneisen  equation  of  state,  and  theoretical 
studies  on  shock  stability  and  phase  transitions.  The  stresses  in  the 
tests  on  wet  playa  ranged  from  1  to  400  kbar.  The  Hugoniot  of  the  wet 
playa  could  be  synthesized  from  those  of  dry  playa  and  water.  There 
appeared  to  be  a  phase  transition  to  stishovite  at  the  high  pressures. 
Reversible  and  irreversible  equations  of  state  for  materials  undergoing 
phase  transitions  were  constructed.  These  were  applied  to  the  transition 
between  quartz  and  stishovite.  The  equations  of  state  were  used  in  wave- 
propagation  calculations.  For  computations,  it  was  found  necessary  to 
employ  artificial  viscosity  during  compression  and  rarefaction  both  in 
order  to  properly  damp  oscillations. 

Petersen,  et  al.  (Ref.  7),  reported  Hugoniot  and  release  adiabat  data 
on  natural  and  reconstituted  tuff,  a  porous,  partially  cemented  mixture 
of  quartz,  feldspar,  pumice,  and  other  minerals.  No  precursor  was  evident 
and  consolidation  appeared  to  occur  at  about  1 00  kbar.  The  release 
adiabats  (unloading  curves)  showed  recovery  to  specific  volumes  larger 
than  the  volume  of  the  solid  material  but  less  than  the  initial  volume 
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even  for  stresses  well  above  100  kbar. 
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b.  Metals 

Kormer  et  al .  (Ref.  8)  conducted  a  series  of  shock  compression  tests 
on  porous  aluminum,  copper,  lead,  and  nickel  in  the  megabar  pressure  range. 
Initial  densities  were  as  low  as  one-fourth  the  crystal  density.  The  work 
was  done  with  bonded  metal  powders  rather  than  sintered  foams.  Their  test 
data  led  them  to  more  complex  formulations  of  the  Mie-Griineisen  equation 
and  to  determination  of  the  variations  of  the  Gruneisen  ratio.  The  Kormer 
equation  has  the  form  of  two  simultaneous  equations,  one  for  pressure  and 
one  for  energy,  both  written  as  functions  of  temperature. 

Krupnikov  et  al .  (Ref.  9)  conducted  shock  compression  tests  for 
porous  tungsten  (initial  densities  as  low  as  one-fourth  the  crystal  density) 
in  the  megabar  range.  They  were  interested  in  formulating  a  refined 

equation  of  state  that  would  include  electronic  components.  The  form 
of  the  equation  is  somewhat  simpler  than  that  of  Kormer  et  al.  and 
could  be  easily  programmed  for  computation.  For  tungsten,  it  was 
found  that  the  effective  value  of  Gruneisen 's  ratio  varied  from  1.5  at 
room  conditions  to  around  .17  at  3.5  Mbar. 

Boade  (Ref.  10)  reported  the  results  of  15  experiments  with  porous 
tungsten  (initially  at  19  percent  of  solid  density)  up  to  pressures  of 
500  kbar.  His  data  gave  a  good  fit  to  the  following  equation: 


U  =  0.0204  +  1.116  U  (6) 

s 

where  U  and  U  are  shock  and  particle  velocity,  respectively,  and  the 
s 

units  are  in  cm/^sec.  The  corresponding  pressure/particle  velocity 
(P  -  U)  plot  also  gave  a  smooth  fit  to  the  data.  However,  the  pressure- 
volume  plot  showed  considerable  scatter. 
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Rempel  et  al .  (Ref.  1)  conducted  shock  attenuation  experiments  on 
several  solid  and  distended  materials.  Tests  on  porous  aluminum 

3  3  3 

(0.7  to  1.4  g/cm  ),  polyuiethane  (.67  g/cm  ),  bervllium  (1.1  g/cm  ) , 
graphite  (1.1  and  1.7  g/cm  )  and  silica  (1.0  g/cm  )  were  mode  in  a 
pressure  range  from  0.1  to  6  kbar.  The  peak  stresses  and  wave  shapes 
were  not  predictable  on  the  basis  of  theory  available  at  the  time.  The 
fore-runner  or  precursor  wave  traveled  at  a  velocity  lower  than  the 

measured  longitudinal  sound  speed  but  higher  than  that  computed  from  the 
measured  elastic  moduli.  The  stress  was  usually  about  uouble  that 
predicted  from  static  compression  tests.  Some  dynamic  and  static 
experiments  were  made  with  preheated  foams.  Preheating  reduced  the 
elastic  modulus,  the  elastic  yield  limit  and  the  stress  level  of  the 
precursor,  and  it  reduced  the  rate  of  attenuation  of  stress  in  the  main 
wave. 

Linde  and  Schmidt  (Ref.  4)  reported  studies  on  porous  aluminum 
(initially  1.08  to  2.14  g/cm3)  and  graphite  (0.9  to  1.7  g/cm3).  The 
experimentally  determined  stress-time  profiles  were  compared  with  wave - 
propagation  calculations  based  on  an  elastic/perfectly  plastic  locking 
model  (shown  in  Fig.  Id).  It  was  found  necessary  to  construct  a  some¬ 
what  more  complicated  Hugoniot  model,  resulting  in  the  AFOAM  computer 
program  for  wave  propagation  calculations.  Comparisons  between  the 
experiments  and  calculations  for  attenuation  tests  showed  that  a 
reasonably  accurate  procedure  for  predicting  attenuation  of  stress  had 
been  developed.  In  some  cases  the  entire  stress  history  could  be 
correctly  predicted. 

Johnson  a,  i  Wackerle  (Ref.  11)  conducted  a  preliminary  study  of 
shock  compression  in  porous  magnesium  and  ammonium  sulphate  as  part  of 
a  study  of  granular  explosives.  The  samples  were  formed  by  pressing 
magnesium  powder  to  75  percent  of  solid  density  and  sulphate  crystals  to 
74  percent  and  98  percent  of  solid  density.  Attempts  to  correlate  the 
experimental  results  with  Thouvenin's  model  (Ref.  34)  were  not  successful, 
but  the  authors  thought  that  the  lack  of  correlation  might  have  been  the 
result  of  experimental  difficulties. 
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Anderson  and  Fahrenbruch  (Ref.  12)  conducted  a  series  of  shock  wave 

| 

experiments  on  porous  and  solid  samples  of  aluminum  and  Teflon  (E.  I.  du 
Pont  de  Nemours  and  Co.).  Aluminum  samples  with  densities  of  1.35,  1.60, 
and  1.92  g/cm3  were  produced  by  cold  pressing  atomized  aluminum  powder. 

The  aluminum  data  from  the  Hugoniot  experiments,  which  were  conducted  with 
shock  stresses  of  200  to  700  kbar,  where  fitted  to  the  equation. 

S  =  £  K  +  A(V)  (7) 


where 

S  =  stress 

T  =  Gruneisen's  ratio 
V  =  specific  volume 
E  =  specific  internal  energy 
A  =  function  of  volume 


The  ratio  f/V,  which  was  found  to  be  a  constant,  was  equal  to  5.15  g/cm3. 
An  expression  for  A  was  also  given. 


Boade  (Ref.  13)  studied  shock  compression  of  porous  copper  (initial 
densities  of  68  and  83  percent  of  solid)  at  low  stresses  and  reported  the 
observation  of  a  three-wave  shock  profile.  The  first  wave  had  an 
amplitude  of  about  one-half  kbar  and  traveled  at  sonic  velocity;  the 
second  had  an  amplitude  of  about  1.1  to  1.3  kbar  traveling  at  about 
one-half  sonic  velocity;  and  the  third  wave  was  the  main  wave.  Complete 
compaction  of  the  material  appeared  to  occur  at  about  21  kbar.  Microscopic 
examination  of  the  materials  showed  that  the  pores  were  essentially 
spherical  and  uniformly  distributed  with  diameters  of  about  0.1  mm.  For 

M  . 

the  stress  levels  used  it  was  found  that  the  ratio  of  Gruneisen  s  ratio 
to  specific  volume  was  essentially  a  constant.  A  Gruneisen  ratio  of 
1.96  was  found  at  the  low  stress  levels.  The  Hugoniot  relation  for  the 
porous  material  was  compared  with  the  formulation  of  Herrmann  (Ref.  44), 
but  the  fit  was  not  very  good.  An  attempt  was  made  to  explain  why  a 
three-wave  structure  might  exist. 
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Butcher  and  Karnes  (Ref.  14)  studied  shock  compression  samples  of 
porous  iron  with  initial  densities  .anging  from  1.3  to  7.0  g/cm3.  Their 
materials  contained  essentially  spherical  voids  0.1  mm  in  diameter  for 
the  lowest  density  foams  and  0.02  mm  for  the  higher  density  sintered 
foams.  The  data  were  interpreted  on  the  basis  of  steady  state  conditions. 
For  the  iron  the  relation  between  the  yield  strength  in  kilobars  and  the 
density  was  found  to  be 


Y  =  8.5  (P  /p  ) 
s  o 


where 

Pg  =  the  density  of  solid  material  at  zero  stress 
=  the  initial  density  of  the  foam. 

The  relation  between  the  initial  wave  velocity  and  density  is  similar 
to  that  between  yield  and  initial  density.  The  compaction  appeared  to 
be  complete  at  26  kbar.  For  correlation  with  experimental  work, 
calculations  were  made  using  either  an  elastic/perfectly  plastic  locking 
model  or  a  simple  locking  model  (see  Fig.  1).  The  simple  locking  model 
was  found  sufficient  for  correlation  with  experiments  on  the  less  dense 
foams.  The  elastic/perfectly  plastic  locking  model  was  necessary  for 
correlation  with  the  higher  density  foams  that  exhibited  a  precursor. 

Schmidt  and  Linde  (Ref.  5)  conducted  a  series  of  tests  on  distended 
copper,  iron,  and  tungsten,  all  at  about  70  percent  of  crystal  density. 
Hugoniots  and  release  states  for  pressures  up  to  about  60  kbar  in  copper, 
50  kbar  in  iron,  and  140  kbar  in  tungsten  were  studied.  It  was  found  that 
inclusions  such  as  corrosion  products  in  the  pores  of  the  porous  materials 
substantially  altered  the  compaction  behavior  of  the  porous  copper  and 
iron.  Clean  samples  of  the  porous  copper  and  iron  compacted  essentially 
to  solid  at  pressures  of  about  20  kbar.  The  tungsten  foam  had  not  fully 
compacted  to  the  solid  state  at  140  kbar. 

c .  Porous  Carbon  Material 

Butcher  (Ref.  15)  reviewed  experiments  on  three  types  of  porous 
carbon  materials  with  initial  densities  of  0.68,  1.14  and  1.30  g/cm3. 
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The  shock  wave  data  on  porous  carbon  materials  were  compared  with  the 
Hugonlot  of  the  solid  for  tests  in  which  the  porous  material  reached 
consolidation  pressures.  Because  of  the  various  types  of  carbon,  and 
because  of  the  phase  transformations  that  occur,  even  the  choice  of  the 
correct  Hugonlot  for  the  solid  was  somewhat  uncertain.  Particle  velocities 
and  wave  velocities  were  compared  with  calculations  based  on  the  plate 
gap  model  of  Thouvenin  (Ref.  35).  The  correlation  between  experimental 
data  and  theoretical  calculations  was  good,  especially  at  the  point  of 
complete  consolidation.  Results  on  the  carbon  foams  show  significant 
attenuation  of  peak  stress  and,  hence,  inapplicability  of  the  usual 
shock-front  relations, which  require  the  presence  of  a  steady  state  in 
the  wave  front . 

Rempel  et  al.  (Ref.  1)  and  Linde  and  Schmidt  (Ref.  4)  conducted 
combined  theoretical  and  experimental  studies  of  graphites:  their 
general  results  have  been  summarized  earlier  under  subsection  b. 

Boade  (Ref.  16)  studied  the  shock  compression  of  foamed  graphite 
with  stresses  in  the  range  from  8  to  190  kbar.  The  initial  densities 
were  0.55  and  0.68  g/cm3.  Hugoniots  were  developed  for  foamed  graphite 
for  several  initial  porosities.  The  data  reduction  was  done  on  the 
assumption  of  a  constant  wave  form  with  a  single  wave  structure. 

However,  in  some  cases  a  two-wave  structure  was  present.  The  precursor 
amplitude  appeared  to  increase  with  the  applied  stress  level. 

d .  Miscellaneous  Materials 

Polyurethane,  Beryllium,  and  silica  were  among  the  materials  studied 
by  Rempel  et  al.  (Ref.  1).  Their  results  were  discussed  under  subsection 

b. 

Johnson  and  Wackerle  (Ref.  11)  made  a  preliminary  study  of  ammonium 
sulphate,  as  mentioned  previously  under  subsection  b. 

Leibermann  (Ref.  18)  conducted  a  series  of  tests  up  to  about  1  kbar 
on  Celotex,  styrofoam,  redwood,  sugar  pine,  foam  glass,  and  balsa.  The 
Hugoniots  obtained  were  similar  to  the  isothermal  equation  of  state  and 
the  data  was  reduced  by  means  of  simple  characteristic  plots. 


'I**.'**.*:  .'"Hywrw re ^»>*W *1* »«»«'»»  W «o»  T>"»  **  *  ,  * 


■  t, 


•mw' 


Lee  (Ref.  19)  studied  dynamic  compaction  of  distended  boron  nitride 
with  an  initial  density  of  about  half  the  solid  density.  The  grain  size 

Q 

was  between  20  and  100  A.  With  pressures  up  to  10  kbar  the  precursor 
amplitude  was  1.5  kbar.  Compaction  was  not  complete  at  the  10  kbar 
pressure  level.  The  experimental  data  were  compared  with  theoretical 
calculations  based  on  the  simple  elastic/perfectly  plastic  locking 
model  and  on  Herrmann's  model  (Ref.  41).  Good  correlation  with  the 
first  model  was  obtained  in  the  P-U  plane  with  the  assumption  of  a 
yield  strength  of  2  kbar. 

Linde  and  DeCarli  (Ref.  20)  studied  phase  transformations  and 
Hugoniot  elastic  limits  in  porous  and  solid  titanium  dioxide  under 
shock  conditions.  For  the  solid,  the  elastic  limit  was  70  or  100  kbar, 
depending  on  the  orientation  of  the  crystal  lattice  in  the  shock  wave. 
Porous  samples  were  formed  to  densities  of  98,  96,  66  and  50  percent  of 
sol*d  density  by  hot  pressing.  Hugoniot  elastic  limits  corresponding  to 
these  densities — 55,  33,  18,  5  kbar — show  a  marked  sensitivity  to 
porosity. 
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SECTION  III 


THEORETICAL  MODEL  FOR  A  POROUS  MATERIAL 
1 .  INTRODUCTION 

The  major  theoretical  contributions  of  this  report  are: 

#  A  set  of  constitutive  relations  for  porous  materials. 

These  relations  provide  for  compaction  to  solid,  and 
for  loading,  unloading,  yielding,  heating,  melting, 
cooling,  and  spalling,  both  before  and  after  full 
compaction  of  the  material  has  occurred. 

•  A  modified  PUFF-type  computer  program  for  calculating  the 
propagation  of  stress  waves  through  a  series  of  solid  or 
porous  materials.  This  modified  code  is  based  on  the  Lax- 
Wendroff  procedure  for  integrating  the  equations  of  motion. 

The  SRI  PUFF  1  code  is  a  computer  program  developed  for  calculating 
one-dimensional  stress-wave  propagation  through  solid  or  porous  materials. 
The  stress  waves  being  computed  are  initiated  either  by  the  deposition 
of  radiated  energy  in  the  materials  or  by  the  impact  of  one  material  on 
another.  Computations  are  made  with  the  Lagrangian  form  of  the  differ¬ 
ential  equations  of  motion,  so  that  coordinates  move  with  the  materials. 
The  program  is  termed  a  Q-code  because  of  the  computation  of  an  artificial 
viscosity  stress,  or  Q,  to  provide  for  the  propagation  of  shock  waves. 


The  following  features  are  included  in  the  code: 

•  Interfaces  between  different  materials  are  allowed;  hence, 
the  stress  wave  propagation  may  be  followed  through  a  series 

of  materials. 

•  The  equation  of  state  of  a  solid  material  is  represented  by 
two  algebraic  equations:  the  Mie-Gruneisen  equation  and  the 
expansion  equation  employed  in  previous  PUFF  versions. 

•  Deviator  stress,  yield  strength,  and  work  hardening  are 

provided  for. 

•  The  materials  may  be  distended  (porous).  Strength  of 
the  porous  material  depends  on  the  internal  energy; 
hence,  it  is  a  function  of  temperature. 

•  Materials  may  spall  or  separate  at  any  coordinate  and 
may  subsequently  recombine.  For  porous  materials,  the 
spall  strength  depends  on  the  degree  of  compaction  and 

on  internal  energy. 
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•  The  radiant  energy  input  is  deposited  over  a  period  of 
time  so  that  material  may  expand  during  deposition. 

ft  The  radiation  absorption  spectrum  for  a  compound 
material  may  be  introduced  as  a  seriis  of  spectra, 
one  for  each  component  or  element  of  the  material. 


A  complete  thermodynamic  equation  of  state  for  a  substance  is  the 
locus  of  all  possible  equilibrium  states  for  that  substance.  Each  state 
is  a  set  of  values  of  the  five  thermodynamic  quantities ;  tensorial  stress, 
specific  volume,  specific  entropy,  specific  internal  energy,  and  temperature. 
Very  shortly  after  a  shock  passes  through  a  homogeneous  solid,  the  material 
usually  reaches  thermodynamic  equilibrium  and  the  macroscopic  thermodynamic 
quantities  vary  smoothly  in  *he  material.  However,  the  states  reached 

within  a  porous  material  are  more  complex.  Stresses  initially  will  be 
higher  at  particle  contact  points,  and  large  strains  and  distortions 
will  occur  as  the  voids  are  filled  in.  Processes  such  as  microjetting 
and  the  production  of  localized  "hot  spots"  v’ill  occur.  For  the  scale 
of  homogeneity  of  most  foams,  thermodynamic  equilibrium  generally  will 
not  be  achieved  within  the  time  scale  of  laboratory  experiments.  Also, 
the  state  reached  by  a  shock  is  not  a  unique  function  of  the  thermo¬ 
dynamic  quantities,  but  depends  on  the  history  of  the  specimen. 


For  the  above  two  reasons  (lack  of  equilibrium  and  lack  of  uniqueness) 
no  true  equation  of  state  can  be  applied  to  the  porous  material  as  a 
whole.  However,  we  can,  and  have,  developed  constitutive  relations  that 
will  serve  the  same  function  as  an  equation  of  state.  In  the  remainder 
of  the  report  we  will  employ  the  popular  terminology  and  loosely  refer  to 
the  constitutive  relations  as  porous  equations  of  state."  In  the  present 
formulation  we  have  accounted  for  the  history  of  the  material.  We  will 
also  assume  that  average  thermodynamic  quantities  are  equilibrium  values 
and  that  when  the  material  is  fully  compacted  to  a  solid,  the  average 
quantities  lie  on  the  true  equation-of-state  surface  of  the  solid. 
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For  convenience  in  formulating  the  equations  of  state  for  both 
solid  and  porous  materials,  the  stress  tensor  will  be  defined  a9  the 
3um  of  a  pressure  and  a  deviator  stress  tensor.  The  pressure  is 
defined  as 

P  =  -  (S  +  S  +  S  )  10) 

J  11  22  33 

where  the  S's  are  stresses  on  any  three  mutually  orthogonal  planes.  The 
deviator  stress  is  the  variation  of  any  normal  stress  from  the  average: 


SD  =  S 
ii  ii 


(i  =  1,  2,  3) 


For  the  one-dimensional  strain  case  to  be  treated  here,  the  stresses  can 
be  conveniently  taken  in  the  directions  of  principal  stresses  so  that  only 
normal  stresses  occur.  Furthermore,  the  deviator  stresses  are  simply 
related, 

SP  =  SD  =  -  \  SD 
aa  33  i  11 

so  that  the  only  stress  quantities  to  be  computed  are  SD  and  P.  The 

ii 

pressure  is  computed  as  a  function  of  two  or  more  of  the  other  thermo¬ 
dynamic  quantities.  The  deviatoric  stress  is  computed  from  a  stress- 

strain  relation.  The  stress  is  then  found  as  a  simple  sum  of  SD  and  P. 

1 1 

With  this  separation  of  stress  into  two  components,  the  development  of  an 
equation  of  state  requires  the  construction  of  two  relationships,  one  for 
pressure  and  one  for  deviator  stress.  This  component  approach  has  been 
U3ed  ;n  the  earlier  development  of  PUFF  and  is  continued  here  for  building 
the  model  for  distended  materials. 

2.  EQUATION  OF  STATE  OF  A  SOLID 

We  consider  first  an  equation  of  state  of  a  solid  (i.e.,  nonporous) 
material.  The  substance  is  termed  solid  because  it  is  solid  at  the 
initial  pressure  and  internal  energy;  however,  the  equation  of  state 
describes  the  material  behavior  in  solid,  liquid,  and  gaseous  phases. 

The  equation  used  here  for  pressure  is  of  the  form 

P  -  P  (E ,  V )  (12) 


I. 
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which  says  that  prossure  is  a  function  of  only  specific  internal  energy 
and  spocific  volume.  The  thermodynamic  quantities  entropy  and  temperature 
are  not  considered  explicitly.  Equation  12  defines  a  surface  in  E-P-V 
space . 

An  equation  of  state  represents  equilibrium  states.  Therefore,  as 
a  material  undergoes  gradual  changes,  such  as  .eating,  compression,  etc., 
the  successive  states  describe  a  path  on  the  equation-of-state  surface  if 
there  is  no  heat  conduction  or  other  nonequilibrium  process  occurring. 

If  the  material  is  compressed  by  passing  through  a  steady-state  shock 
front  and  the  initial  and  final  states  are  equilibrium  states,  then 
these  states  lie  on  the  equation-of-state  surface.  These  initial  and 
final  states  are  connected  by  a  straight  line,  the  Rayleigh  line,  which 
lies  on  or  above  the  surface,  for  the  usual,  concave-upward,  surfaces. 

The  states  of  transition  within  a  shock  front  are  not  states  of  thermo¬ 
dynamic  equilibrium  and  hence  do  not  necessarily  lie  on  the  surface. 

Shock  experiments  lead  to  the  determination  of  a  Hugoniot  or  Rankine- 
Hugoniot  equation  of  ata-e  which  i3  not  a  complete  equation  of  state,  but 
represents  one  curve  on  the  equat j on-of-state  surface.  This  line  is  the 
locus  of  final  states  that  can  be  obtained  by  a  steady-state  shock 
transition  from  a  given  initial  state.  A  common  form  for  a  pressure-volume 
Fugonirt  is  shown  in  Fig.  3, 

This  same  curve  is  redrawn  in  Fig.  4  on  the  E-P-V  equation-of-state 
surface  for  the  material.  During  compression  there  is  some  increase  in 
internal  energy,  so  that  the  Hugoniot  does  not  lie  in  a  single  P-V  plane. 

As  a  reminder  of  the  role  of  stress  in  the  compression  of  the  solid, 
consider  the  stress-volume  Hugoniot  of  Fig.  3.  During  compression  the 
stress  is  greater  than  the  pressure;  on  unloading,  the  stress  decreases 
rapidly  to  yielding  and  then  follows  a  stress  adiabat  below  the  pressure 
adiabat.  The  unloading  (pressure)  adiabat  is  also  depicted  in  Fig.  4. 

The  unloading  adiabat  lies  to  the  right  of  the  Hugoniot  for  materials 
that  expand  during  heating.  For  such  materials  less  internal  energy 
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is  recoverod  by  tho  adiabatic  decompression  than  was  insorted  to 
produce  the  shock  compression.  A  line  illustrating  a  path  of  adiabatic, 
gradual  compression  is  also  shown  in  Fig.  4. 

Several  other  lines  of  interest  are  shown  in  Fig.  4.  The  zero 
(or  atmospheric)  pressure  line  is  the  locus  of  points  obtained  by  simply 
heating  the  material  without  external  mechanical  confinement.  Heating 
increases  the  internal  energy,  and  thermal  expansion  occurs.  For  small 
increases  in  internal  energy,  the  zero  pressure  curve  describes  the  usual 
expression  for  volumetric  thermal  expansion 

V  =  V0  (1  +  OAQ  ) 

or 


where 

Ot  =  the  volumetric  thermal  expansion  coefficient 
f  =  the  Gruneisen  ratio 
K  =  the  bulk  modulus 
A0  =  the  change  in  temperature 
AE  =  the  change  in  internal  energy 
V0  =  the  initial  specific  volume 

If  we  take  the  usual  assumption  that  F/V  and  K  are  constant  at  zero 
pressure,  then  AE  is  proportional  to  A 0  while  the  material  is  solid. 

As  the  material  is  further  heated  it  melts  and  then  vaporizes.  The 
P  =  0  curve  becomes  asymptotic  to  the  line  described  by 

E  =  sublimation  energy 
P  =  0 

for  large  V. 

The  region  of  spalling  is  not  well-defined  at  present,  so  that 
the  depicted  surface  below  zero  pressure  is  conjectural.  Spalling,  or 
mechanical  separation  of  the  material,  occurs  if  the  solid  material  is 
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subjected  to  a  suffici  ntly  large  tensile  stress.  The  occurrence  of 
spalling  depends  on  the  magnitude  of  the  applied  stress  and,  in 
general,  on  all  three  principal  stresses  at  a  point.  (Hence,  there 
is  considerable  ambiguity  in  depicting  spalling  on  a  pressure  coordinate 
as  in  Fig.  4.)  Because  of  the  time  required  to  form  and  propagate 
cracks  and  to  move  masses  of  material  so  that  separation  can  occur,  the 
duration  of  loading  is  also  significant.  The  spall  resistance  will  also 
depend  on  temperature  or  internal  energy.  If  the  material  is  liquid  at 
the  time  of  separation  and  the  separation  is  accompanied  by  the  production 
of  bubbles  rather  than  the  growth  of  cracks,  the  phenomenon  is  termed 
cavitation.  Bull  (Ref.  50)  has  pointed  out  that  the  occurrence  of 
cavitation  is  governed  by  the  size  of  cavitation  nuclei  (bubbles  or 
foreign  particles),  surface  tension,  viscosity  and  inertia  forces,  and 
the  amplitude  and  duration  of  the  applied  stress.  Of  these,  the  viscosity 
and  stress  amplitude  will  probably  dominate  under  shock  loading.  When 
the  material  reaches  the  point  of  vaporization, it  ceases  to  have  any 
resistance  to  cavitation.  For  the  SRI  PUFF  1  code  a  single  tensile  stress 
value  is  used  for  the  spallation  criterion  at  present. 

When  spelling  occurs,  the  stress  on  the  spall  surface  immediately 
returns  to  zero.  The  material  some  distance  from  the  spall  surface  is 
relieved  gradually  by  rarefaction  waves  from  the  spall  surface.  This 
deeper  material  then  returns  to  higher  compression  stresses  upward 
along  a  path  similar  to  the  original  spall  path,  recovering  internal 
energy  as  it  recompresses.  The  region  to  the  right  of  the  spall 
criterion  cannot  be  reached.  It  may  be  noted  that  specific  volume  is 
taken  to  mean  the  total  volume  of  solid  particles  on  both  sides  of 
the  spall  (but  not  including  the  volume  of  the  space  between  the  spall 
surfaces)  divided  by  the  mass  of  the  particles. 

When  material  is  held  at  a  particular  volume  and  heated  (internal 
energy  is  added),  it  goes  through  states  that  are  straight  lines  on  the 
equation-of -state  surface.  This  is  indicative  of  the  fact  that  for 
constant  V (=  V  )  the  analytical  equations  for  the  surface  have  the 

l 

form 
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E  =  A(V  )  •  PV 
1  1 


where  A(V  )  is  a  function  of  V  only. 
1  l 


The  equation-of-state  surface  depicted  in  Fig.  4  is  an  idealized 
form  applicable  to  a  material  that  does  not  experience  phase  changes 
or  other  phenomena  that  lead  to  regions  of  negative  curvature  in  the 
P-V  plane.  While  this  surface  represents  the  material  behavior  quali¬ 
tatively,  only  certain  regions  of  the  surface  are  well  understood 
quantitatively.  The  best  understood  region  is  in  the  vicinity  of  the 
Hugoniot  because  of  the  availability  of  experimental  data  along  that 
curve.  The  least  understood  regions  are  those  near  spalling  and  Just 
to  the  right  of  the  curve  V  *  VQ. 

The  equation  of  state  for  a  solid  and  the  subroutine  for  calculating 
it  in  SRI  PUFF  1  are  essentially  the  same  as  those  in  PUFF  66.  The 
equation  of  state  is  described  by  two  analytical  equations  (one  for 
compression  and  one  for  expansion)  and  is  bounded  for  negative 
stresses  by  a  spall  criterion. 

The  equation  used  to  describe  compression  is  the  Mie-Gruneisen 


equation 
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where 


P  and  E  =  a  point  on  some  reference  curve  at  the  same 
REF  REF 

specific  volume  V 
r(V)  =  the  Gruneisen  ratio. 

Equation  (15)  has  been  derived  on  the  assumption  that  T  is  a 
function  of  V  only.  Equation  (15)  provides  a  means  for  extending  the 
information  of  a  known  P-V  relation  (such  as  a  Hugoniot)  to  other  values 
of  internal  energy.  Because  the  Hugoniot  is  the  P-V  relation  that  is 
most  likely  to  be  known,  the  computations  are  constructed  so  that  the 
Hugoniot  is  the  reference  curve  used.  The  Hugoniot  P-V  equation  is 
presumed  to  be  in  the  form 


*<*<**■ 


P„  =  Cp  +  Dp2  +  S|i3  (16) 

H 

where 


The  internal  energy  along  the  Hugoniot  is 

EH  =  *PH<V°  -  VH>  (17) 

Equation  (17)  is  based  on  the  assumption  that  the  initial  internal 
energy  is  zero  and  that  the  Hugoniot  is  everywhere  concave  upward.  In 
general,  tne  latter  assumption  excludes  consideration  of  phase  changes. 
Also,  the  relation  is  strictly  true  only  if  the  stress  Hugoniot  coincides 
with  the  pressure  Hugoniot;  however,  at  high  pressures  there  is  usually 
little  inaccuracy  introduced  by  this  approximation.  With  the  aid  of 
Equations  (16)  and  (17)  the  Mie-Gruneisen  equation  takes  the  following 
form  in  the  program:  ,  . 

P  =  (Cp  +  Djl2  +  Sp3)  *  (l  -  TpE  (18) 

In  the  computer  program  the  Gruneisen  ratio  T0  at  initial  density 
is  taken  as  a  constant,  EQSTG.  Then  F  is  treated  as  a  function  of 
density  such  that  Tp  is  constant. 

At  a  constant  volume,  Equation  (18)  has  the  form  of  Equation  (14); 
hence,  constant  volume  lines  on  the  equation-of -state  surface  are 
straight  lines.  The  Mie-Gruneisen  equation  of  state  is  used  for 
densities  greater  than  the  initial  density.  Thus  on  the  equation-of - 
state  surface  the  straight  line  V  =  V0  is  the  boundary  between  the 
Mie-Griineisen  equation  and  an  expansion  equation. 

The  expansion  equation,  which  is  unchanged  from  PUFF  66,  meets  four 
requirements: 

•  It  joins  smoothly  to  the  Mie-Gruneisen  equation  along  V  =  V0. 

•  It  expands  like  PV  =  (y  -  1)E  at  large  expansions  (like  a  gas). 

•  It  provides  a  linear  relation  between  P  and  E  for  constant  V. 
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•  It  accounts  for  the  partition  of  internal  energy  into 
components  for  kinetic  energy  and  for  intermolecular 
bond  disintegration  (sublimation). 

An  equation  which  satisfies  these  requirements  is 


P  = 


where 


H  +  (T  -  H)  — 
Po 


*  *  |e  -  Eg  1  -  exp(N  (1-Po^  )p0^ 


p  =  density 
p0  =  initial  density 
F  ~  Gruneisen  ratio 

H  =  v  -  1  =  C„/C  -  1  for  expansion  at  low  densities 

'  P  v 


N  =  C/  (TEgPo ) 


S  =  sublimation  energy 

C  =  coefficient  in  Equation  (18),  the  bulk  modulus  at 
low  pressures 

In  the  PUFF  66  manual  (Ref.  35),  a  value  of  0.25  is  suggested  for  H. 
The  sublimation  energy  as  defined  there  is  the  difference  between  the 
internal  energy  of  the  solid  material  at  ambient  conditions  and  the 
internal  energy  of  the  fully  expanded  vapor  at  a  temperature  of  absolute 


The  form  of  the  equation-of -state  surface  generated  by  the  two 
equations  and  the  spall  criterion  is  shown  in  Fig.  4.  When  the  spall 
criterion  is  omitted,  cuts  through  the  surface  at  constant  energy 
produce  the  curves  shown  in  Fig.  5  (for  aluminum).  For  internal  energy 
less  than  the  sublimation  energy,  these  curves  indicate  that  as  the 
material  is  expanded  the  pressure  decreases  to  some  minimum  and  then 
increases  again  to  zero.  In  many  cases  this  uniform  expansion  of 
continuous  material  will  be  interrupted  by  spalling  or  cavitation. 
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If  the  material  is  heated  without  external  confinement,  the  path 
is  on  the  intersection  of  the  P  =  0  surface  and  the  equation-of -state 
surface  (see  Fig.  6).  The  initial  expansion  follows  the  usual  law  of 
thermal  expansion,  Equation  (13). 

The  deviator  stress  equation  takes  a  much  simpler  form  in  the 
PUFF  formulation  than  does  the  pressure  equation.  The  deviatorlc  stress 


Po 


for  SD  <  —  Y 

J 


otherwise 


where 


-  i* 


“  the  shear  modulus 


=  the  density 
=  the  yield  strength 
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FIG  6  RELATION  BETWEEN  ENERGY  AND  VOLUME  FOR  UNCONFINED 
EXPANSION  OF  SOLID 


32 


p-.ww  ffi.TiiTtwiiv *  r 


3.  CONSTITUTIVE  RELATIONS  FOR  A  DISTENDED  MATERIAL 

a .  Introduction 

The  model  proposed  here  for  the  description  of  the  state  of  a 
distended  materia]  consists  of  a  set  of  constitutive  relations  to  define 
behavior  while  the  material  is  porous,  plus  an  equation  of  state  for 
the  solid  substance  to  be  used  in  describing  behavior  after  consolidation. 
As  discussed  above,  the  model  for  behavior  of  porous  material  is  not, 
strictly  speaking,  an  equation  of  state  because  not  only  is  thermodynamic 
equilibrium  generally  not  achieved,  but  there  is  no  unique  relation 
among  E,  P,  and  V.  Instead,  the  pressure  depends  uot  only  on  E  and  V, 
but  also  on  the  previous  E-P-V  states. 

A  Hugoniot  for  a  porous  material  is  shown  in  Fig.  7.  This  Hugoniot 
plus  a  Hugoniot  for  the  same  material  starting  at  the  normal  solid 
density,  are  traced  on  the  equation-of-state  surface  of  Fig.  8.  A 
comparison  of  these  two  Hugoniots  illustrates  that  for  the  same  pressure 
or  volume,  points  on  the  Hugoniot  for  the  porous  material  are  at  higher 
energy  states  because  a  greater  amount  of  compressional  work  (J'pdV)  is 
expended  to  reach  a  given  pressure  or  volume.  The  individual  rounded 
particles  undergoing  shock  compression  are  loaded  by  a  highly  nonuniform 
stress  field,  are  caused  by  yield  and  flow  to  fill  the  voids,  and  are 
locally  heated  by  the  energy  of  deformation  and  other  processes  (e.g., 
microjetting  and  adiabatic  compression  of  any  gas  in  the  pores) . 

The  portion  of  the  constitutive  relations  that  describes  the 
behavior  of  the  unconsolidated  material  with  solid  particles  is  the  main 
feature  to  be  developed  (see  Fig.  9).  Any  point  in  this  three-dimensional 
region  represents  a  possible  state  for  the  material. 


* 

This  point  is  discussed  furv.her  in  subsection  c  below. 
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FIG.  7  HUGONIOT  FOR  A  POROUS  MATERIAL 


\ 


NOT  REPRODUCIBLE 


HUGONiOT 


YIELD  SURFACE  OF 
CONSTITUTIVE  RELATIONS 


COMST  f.N  r  VOLUME 


^UGt.Nlv*-  OF  POROUS 
MATLrli  AL 


EQUATION -OF -STATE 
SURFACE  FOR  SOLID 


UNCONFINED 

thermal 

EXPANSION 


OF  SOLID 


•  P-IIH-  t( 


FIG.  8  DEPICTION  IN  E-P-V  SPACE  OF  THE  CONSTITUTIVE  RELATIONS 
FOR  A  POROUS  MATERIAL 
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FIG.  9  CONSTITUTIVE  RELATIONS  OF  A  POROUS  MATERIAL,  EMPHASIZING 
THE  INTERMEDIATE  SURFACE  FOR  ELASTIC  RESPONSE  TO 
LOADING  AND  HEATING 

The  region  is  bounded  above  by  a  surface  prescribed  by  the  yielding  of 
the  particles  (surface  EBGFCH) ,  below  by  spalling  or  separation  of  the 
particles  (surface  EADFCH) ,  and  to  the  left  by  consolidation  to  a  solid 
material  (surface  DFG) .  When  the  state  point  (R,  as  an  example)  is 
within  this  region,  the  state  point  is  constrained  to  move  on  an  inter¬ 
mediate  surface  (ABC)  which  is  approximately  parallel  to  the  equation- 
of -state  surface  for  a  solid  at  the  same  E  and  P  values.  The  detailed 
character  of  the  intermediate  and  bounding  surfaces  is  developed  in  the 
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following  paragraphs,  based  on  a  consideration  of  various  processes  that 
material  may  undergo:  loading,  unloading,  and  reloading:  heating  and 
cooling;  and  spalling. 

The  behavior  of  the  porous  material  under  heating  at  zero  pressure 

can  be  estimated  by  considering  the  behavior  of  the  solid  particles.  If 

the  material  is  heated  from  an  arbitrary  initial  point  at  zero  pressure, 

* 

such  as  point  R  in  Fig.  9,  the  expansion  path  is  along  the  curve  RN 
defined  by 


P  =  0 


v  =  Vj  (1  +  a  L  9) 


where  £0  is  a  change  in  temperature.  Note  that  this  expression  for 
the  path  holds  for  Vx  =  V0  ,  the  initial  specific  volume,  and  V1  =  Vs , 
the  solid  volume.  Because  this  path  represents  no  mechanical  yielding, 
particle  rearrangement,  melting,  or  spalling,  it  is  presumed  to  lie  on 
an  intermediate  surface  ABC  of  Fig.  9. 


Unloading  from  a  partially  consolidated  state  will  occur  along  a 
different  path  than  that  along  which  the  loading  (including  nonshock 
loadings)  has  occurred;  i.e.,  the  load-unload  process  is  irreversible 
(nonisentropic) .  A  possible  loading-unloading  path  is  shown  in  Fig.  7. 
This  unloading  path  serves  to  further  define  the  intermediate  surface 
ABC  of  Fig.  9.  The  unloading  path  is  not  in  a  single  P-V  plane  (see  lines 
R’L1  and  RL  in  Fig.  9),  but  slopes  in  the  E  direction  to  indicate  a 
decrease  in  internal  energy  with  unloading.  The  fact  that  ^n^oading 
does  not  return  the  material  to  its  initial  state  is  anticipated  from 
a  consideration  of  the  deformed  particles  of  the  material.  During 
loading, they  have  been  deformed  to  fill  voids.  On  unloading,  they  will 
exhibit  a  rebound  that  will  usually  be  much  smaller  than  the  initial 


*  For  the  computations ,  the  change  in  temperature  is  assumed  to  be 

proportional  to  the  change  in  internal  energy  while  tht  particles  are 
in  the  solid  phase.  Hence  RN  is  a  straight  line  in  the  computations. 
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loading  deformation.*  Hence  the  much  steeper  slope  of  the  unloading  than 
the  loading  P-V  relation  of  Fig.  7.  Reloading  will  occur  mainly  through 
elastic  deformation  of  the  particles  again  and  will  proceed  along  a  path 
similar  to  the  unloading  curve  until  gross  plastic  flow  begins  when  the 
virgin  loading  curve  Is  reached. 

Shear  and  bulk  moduli  along  the  unloading  and  reloading  paths  are 
assumed  to  be  a  function  of  density.  The  moduli  are  assumed  to  vary  in 
the  same  manner  as  the  pressure  along  the  porous  Hugoniot.  Hence,  normally 
the  moduli  will  be  small  initially,  and  correspond  to  the  precursor 
velocity.  The  moduli  will  then  increase  with  compaction,  reaching  the 
value  of  the  solid  moduli  at  complete  consolidation.  The  unloading 
moduli  govern  the  arrival  time  of  small  reflected  waves  that  reach  the 
gage  later  than  the  precursor  but  earlier  than  the  main  wave.  Hence,  by 
comparing  experimental  records  with  computed  histories  it  is  possible  to 
determine  the  variation  of  the  moduli. 

The  behavior  of  the  porous  material  under  combined  loading  and 
heating  determines  the  upper-bound  surface  on  the  constitutive 

relations  of  Fig.  9.  Assume  that  by  some  combination  of  heating  and 
loading  the  state  is  some  point  on  the  P  =  0  surface,  on  the  curve 
V  =  Vj  (1  +  Ot  A  6) ,  Then  let  the  material  be  loaded  adiabatically .  The 
paths  RL  and  R'L'  on  the  intermediate  surface  ABC  of  Fig.  9  represent 
such  loadings.  When  yielding  occurs,  the  particles  flow,  changing  the 
configuration,  and  the  state  point  moves  out  of  the  plane  ABC  along  the 
path  LM.  Yielding  of  the  aggregate  depends  on  the  stresses  and  the 
yield  strengths  of  the  individual  particles.  In  turn,  yield  strength 
depends  on  the  temperature.  Available  data  on  the  variation  of  yield 
strength  with  temperature  are  used  to  construct  the  yield  locus  BC. 

Fig.  10  presents  yield  strengths  of  some  aluminum  alloys  as  a  function 
of  temperature.  The  assumption  is  made  here  that  the  pressure  along 
curve  BC  in  Fig.  9  is  proportional  to  the  thermal  strength  reduction 
curve  of  Fig.  10. 

*  This  is  true  at  early  times  even  for  materials  like  ATJ  graphite, 
which  recover  their  initial  distention  at  some  later  time  (Ref.  4). 
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FIG.  10  VARIATION  OF  STRENGTH  WITH  TEMPERATURE  FOR  1100  ALUMINUM 


If  unloading  proceeds  below  the  P  ■  0  surface,  spalling  may  occur. 
The  strengtn  of  the  interparticle  bonds  will  probably  depend  on  the 
amount  of  precompression  and  heating.  Hence,  spall  strength  should  be  a 
function  of  porosity,  particle  size  and  shape,  degree  to  which  inter¬ 
particle  welding  has  occurred,  as  well  as  on  the  usual  quantities  such 
as  strain  rate,  temperature,  maximum  stress  and  deviatoric  stress,  and 
the  loading  history.  Then  the  spall  strength  forms  a  surface  (or  family 
of  surfaces  if  all  factors  are  considered)  in  the  E-P-V  space.  A 
possible  spall  criterion  surface  is  shown  in  Fig.  9.  Direct  tests  of 
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spall  strength  of  porous  materials  are  not  available.  However, 
estimates  of  such  strength  may  be  made  by  considering  the  extent  to 
which  interparticle  welding  has  occurred  and  by  estimates  from  photo¬ 
micrographs  of  the  areas  of  interparticle  bonds.  The  spall  criterion 
surface  is  defined  for  the  purposes  of  the  computer  program  by  a  linear 
relation  in  the  P-V  plane  from  D  to  E  in  Fig.  9.  The  remainder  of  the 
surface  is  defined  by  its  intersection  with  an  intermediate  surface  ABC. 
The  thermal  strength  reduction  function  of  Fig.  10  is  presumed  to  govern 
tensile  as  well  as  compressive  strength.  With  this  definition  of  the 
spall  surface,  the  behavior  of  the  porous  material  with  particles  in  a 
solid  phase  is  defined. 

The  behavior  of  porous  material  comprised  of  melted  or  gaseous 
particles  will  now  be  examined  by  using  the  surfaces  shown  in  Fig.  11. 
Consider  a  porous  sample  held  at  a  constant  volume  VB  and  initially  at 
E  =  0,  P  =  0.  If  the  material  is  heated  (  in  some  type  of  instantaneous 
flash-heating  process^  the  pressure  increases  as  the  material  tries  to 
expand.  Simultaneously,  the  yield  strength  decreases  with  heating. 

With  further  heating,  yielding  occurs  (point  A  in  Fig.  11),  the 
supported  external  pressure  decreases,  and  the  thermal  expansion  of 
the  particles  partially  fills  the  voids.  When  the  particles  melt 
(point  B) ,  no  external  pressure  is  supported  and  the  molten  particles 
are  in  an  unconsolidated  state  like  water  vapor  in  the  air.  As  heating 
continues , the  particles  expand  and  flow  at  interparticle  contact  points 
until  all  the  voids  are  filled  (point  C)  .  Under  further  heating 

(following  the  filling  of  all  voids)  at  the  constant  volume  V  ,  pressure 

2 

increases  to  some  point  D  in  proportion  to  the  increase  in  internal 
energy  in  accordance  with  the  linear  relation 

E  =  A  (V  )  .  PV  (  I'c) 

2  2 
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FIG.  '1  LOCUS  OF  STATES  CAUSED  BY  HEATING  POROUS 
MATERIAL  WITHOUT  EXPANSION 


Note  that  at  state  point  D,  the  material  is  now  on  the  equation-of -state 


surface  for  the  solid.  The  joint  actually  occurs  over  the  whole  surface 


DFG  and  along  the  P  =  0  line  in  Fig.  9  because  material  may  pass  from 


a  porous  to  a  solid  state  at  any  point  on  the  surface  or  line.  However, 


the  most  commonly  used  consolidation  paths  will  include  points  on  the 


P  =  0  line  and  on  the  line  FG  joining  the  yield  surface  and  the  solid 


equa^lon-of-state  surface.  For  use  in  the  program,  the  joint  line  FG 
is  required  to  pass  through  the  points  at  E  =  0  and  at  E  =  energy  of 


melting.  The  joint  provides  continuity  of  the  two  surfaces  being  joined 


and  does  not  require  a  local  distortion  of  either  surface. 


mmmm 
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b.  Hugoniot  Functions 

There  are  provisions  in  the  code  for  both  P-V  and  SD-V  relations 
for  a  porous  material.  These  relations  are  not  inserted  as  a  single 
function  but  as  a  series  of  functions,  each  pertaining  to  an  interval 
on  the  abscissa  (specific  volume).  Figure  12  exhibits  some  possible 
variations  of  pressure  and  yield  strength.  The  pressure  is  taken  as  a 
series  of  parabolic  segments,  and  the  yield  strength  is  taken  to  vary 
linearly  within  each  specific  volume  interval. 

The  pressure  function  is  not  a  Hugoniot,  as  it  lies  entirely  in 
the  E  -  0  plane.  For  impacts  that  just  produce  consolidation,  the 
internal  energy  is  fairly  small,  so  that  the  Hugoniot  and  the  pressure 
function  will  not  differ  greatly. 

The  yield  strength  information  is  provided  as  a  series  of  linear 
segments  (see  Fig.  12),  each  segment  corresponding  to  one  of  the  porous 
regions  defined  for  the  pressure  Hugoniot.  Both  the  yield  strength  and 
the  shear  modulus  decrease  with  increasing  internal  energy.  Hence,  for 
the  porous  material,  the  deviator  stress  is  also  treated  as  a  function 
of  internal  energy. 
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c  .  Intermediate  Surface  Computations 

The  intermediate  surfaces  are  generated  in  such  a  way  that  loading 
and  unloading,  heating,  and  cooling  are  provided  for  in  a  way  similar  to 
that  for  the  solid  equation-of-state.  The  requisite  E-P-V  equation  is 


P 


Po  " 


K  *  AV 


K  ♦  T  *  AE 
V  *  K 


(23) 


where 

P0  =  the  pressure  at  the  previous  cycle 

V  =  the  specific  volume  at  the  previous  cycle 

K  -  the  b’jl.t  modulus  on  the  intermediate  surface, 

generally  a  function  of  E  and  V 

AV  =  the  change  in  specific  volume 

F  =  the  Gruneisen  ratio 

AE  =  the  change  in  internal  energy  from  the 
previous  cycle 

V  =  the  specific  volume  of  the  solid  particles 
s 

=  the  bulk  modulus  of  the  solid  material 

The  second  term  on  the  right  side  of  Equation  (23)  provides  for 
loading  or  unloading.  The  third  term  provides  for  expansion  or 
contraction  under  heating  and  cooling.  If  P  is  held  constant  in 
Equation  (23)  and  the  material  is  heated,  then 


v  r 

AV  =  — -  AE 

V  Kn 
s  J 


(24) 


For  the  solid  material  at  low  stresses  the  corresponding  relation  is 


AV 


(25) 


Hence,  for  the  same  amount  of  energy  deposition  the  thermal  expansion 

of  the  porous  material  is  larger  than  that  of  the  solid  by  the  factor 

V_/V  ,  as  expected. 
u  s 
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d .  Strength  Redu  t i on 

As  a  material  is  heated.it  loses  strength,  reaching  approximately 
zero  strength  at  melting.  (Curves  showing  this  strength  reduction  are 
shown  in  Fig.  10  for  three  tempers  of  1100  Aluminum.) 

For  computational  purposes,  the  strangth  reduction  function  is 
assumed  to  consist  of  two  parabolas,  as  shown  in  Fig.  13. 


FIG.  13  THERMAL  STRENGTH  REDUCTION  FUNCTIONS  EMPLOYED  IN 
COMPUTER  CODE 


e .  Discussion  of  the  Postulated  Constitutive  Relations 

As  mentioned  above,  it  is  not  strictly  legitimate  to  call  the  postu¬ 
lated  E-P-V  surfaces  for  a  porous  material  an  equation-of-state  and  it 
does  not  appear  that  a  true  equation-of-state  for  the  porous  materials 
is  applicable.  In  the  porous  material,  pressure  is  not  a  unique  function 
of  energy  and  volume,  but  depends  on  the  previous  history  of  the  material. 
The  stresses  and  temperature  vary  throughout  the  porouu  materiel  so  that 
the  constitutive  relations  do  not  represent  equilibrium  states  for  the 
material,  but,  rather,  averages  of  the  thermodynamic  variables.  This 
nonequilibrium  condition  will  persist  for  some  time  after  the  shock  has 
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consolidated  the  material,*  so  that  the  equation  of  sttfte  used  for  the 
solid  will  really  represent  only  effective  average''  states.  Such  averages 
of  the  thermodynamic  variables  do  not  determine  points  that  lie  on  the 
equation-of-state  for  the  solid  particles.  However,  the  average  thermo¬ 
dynamic  quantities  should  lie  very  near  the  surface.  The  strongest  varia¬ 
tion  of  any  variable  within  particles  is  in  E,  the  weakest  variation  is 
in  V.  The  equation-of-state  surface  is  linear  in  E  and  P  for  constant  V, 
and  hence  all  of  the  states  of  the  particles  at  any  instant  lie  near  a 
surface  that  is  approximately  planar.  Thus  a  reasonable  approximation  is 
that  the  average  values  of  the  thermodynamic  variables  lie  near  the  surface. 

Uniqueness  of  the  bounding  surfaces  (yield,  spall,  solidification) 
of  the  constitutive  relations  has  been  assumed  for  the  construction  of 
the  constitutive  relations.  For  the  actual  material  these  surfaces  might 
be  different  if  they  were  reached  by  different  processes  (for  example, 
loading  and  heating,  instead  of  heating  and  loading). 

The  intermediate  surfaces  for  porous  material  (such  as  ABC  in  Fig.  9) 
have  been  somewhat  arbitrarily  constructed  with  a  bulk  modulus  that  varies 
with  temperature  in  the  same  manner  as  the  yield  strength,  and  with  a 
constant  expansion  coefficient.  The  warped  intermediate  surface  thus 
defined  has  not  been  verified  experimentally.  The  path  across  this  surface 
on  loading  (path  LM  in  Fig.  9)  is  like  that  of  an  ideally  plastic  material 
with  work  hardening.  The  actual  path  (and  hence  the  actual  intermediate 
surface)  would  be  rounded  near  yielding  and  might  not  even  contain  a 
linearly  elastic  region.  These  questions  on  the  legitimacy  of  the  E-P-V 
surface  for  a  porous  substance  indicate  that  calculations  made  bv  using 
the  surface  are  not  rigorously  Justifiable  on  thermodynamic  grounds.  The 
surface  may  be  used  for  computations  with  the  realization  of  its  approxi¬ 
mate  nature.  The  equation-of-state  surface  for  a  porous  material  is  some¬ 
what  idealized  but  is  intended  to  account  approximately  for  all  of  the 
important  phenomena.  Because  the  surface  was  generated  from  considerations 
of  temperature  (although  energy  was  used  as  a  variable  rather  than  temper¬ 
ature),  volume,  and  pressure,  those  quantities  are  probably  represented 

The  time  required  for  equilibration  will  be  a  function  of  the  particle 

size. 
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best.  No  attempt  has  been  made  to  indicate  the  variation  of  entropy 
in  the  porous  media. 
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SECTION  IV 


EXPERIMENTAL  MEASUREMENTS 

1.  QUASI-STATIC  MEASUREMENTS  OF  PRESSURE -VOLUME 

Quasi-static  one-dimensional  compression  tests  were  performed  on 
samples  of  the  porous  (70  percent  of  solid  density)  copper,  iron,  and 
tungsten  studied  previously  by  Schmidt  ano  Linde  (Ref.  5).  The  purpose 
of  these  tests  was  to  measure  a  loading  isotherm  and  unloading  and  re¬ 
loading  curves  for  the  three  materials.  It  was  hoped  that  these  measured 
curves  would  then  aid  in  the  formation  of  constitutive  relations  to  des¬ 
cribe  the  behavior  of  the  porous  materials  in  the  part. ally  compacted 
region,''  particularly  just  above  the  elastic'  limit,  where  accurate 
dynamic  measurements  are  difficult  to  obtain.  A  major  question  of 
interest  was  whether  the  quasi-static  P-V  rata  could  be  used  to  approxi¬ 
mate  the  dynamic  data  in  the  partial  compaction'  region,  despite  the 
drastic  differences  in  loading  rate. 

The  testing  apparatus,  shown  schematically  in  Fig.  Id,  was  made  of 
high-strength  steel.  The  capacity  was  limited  to  about  10  kbar.  The 
porous  samples  used  were  about  1  inch  long.  They  were  inserted  in  the 
die  inside  a  1-mil-thick  sleeve  of  indium — the  indium  being  required  for 
lubrication  to  reduce  wall  frictior.  Pressure  was  applied  wi th  a 
hydraulic  Jack  acting  on  the  plunger  above  the  sample.  Tests  were 
performed  by  loading  to  a  certain  pressure  and  reading  the  length  of  the 
plungers  and  sample  with  a  gage  recording  to  0.0001  inch.  Then  the 
pressure  was  released  until  the  gage  reading  began  to  change.  The  "true'' 
pressure  was  then  taken  as  the  average  of  the  upper  and  lower  pressures, 
this  double-reading  procedure  was  employed  to  reduce  errors  associated 
with  friction  between  the  sample  and  the  body  of  the  die. 

To  determine  both  the  loading  isotherm  and  the  unload-reload 
behavior,  two  testing  plans  were  adopted.  In  the  first,  the  loading  was 
increased  monotonically  to  10  kbar.  In  the  second  plan,  the  loading 
proceeded  only  to  some  fraction  of  10  kbar  and  then  the  sample  was 
unloaded  and  reloaded  to  a  higher  pressure.  This  load-unload  cycle  was 
repeated  until  10  kbar  was  reached. 
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FIG.  14  CUTAWAY  VIEW  OF  QUASI-STATIC 
COMPRESSION  DEVICE 


Compression  of  the  sample  was  deduced  from  the  recorded  overall 
shortening  of  the  combination  of  the  sample  and  two  plungers.  In  deter¬ 
mining  the  total  volume  change  of  the  sample,  an  allowance  was  also  made 
for  the  lateral  expansion  of  the  die.  This  lateral  expansion  (  which 
necessitated  a  small  correction  to  the  longitudinal  volume  change)  was 
deduced  fror  calibration  tests  with  sodium  chloride  as  a  sample  with  a 
known  compression  isotherm. 

A  comparison  of  the  results  from  the  first  and  second  test  plans 
indicated  that  the  isotherm  of  the  porous  sample  was  not  altered  signi¬ 
ficantly  by  the  unload-reload  cycle.  Unloading  and  reloading  curves  had 
slopes  similar  to  that  of  nonporous  material;  the  accuracy  of  the  meas¬ 
urements  was  not  sufficient  to  determine  the  variation  of  modulus  with 
density.  The  loading  isotherms  develop  f,  in  three  tests  are  s1  own  in 
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Figs.  15  through  17.  The  irregularities  in  the  curves  shown  are 
associated  with  errors  introduced  in  the  load-ur:load  cycles  and  are 
not  significant  features  of  the  data.  Hugoniot  data  from  Schmidt  and  Linde 
(Ref.  5),  Boade  (Ref.  7,  8)  and  Butcher  and  Karnes  (Ref.  11)  are  also 
shown  in  these  figures.  It  is  noted  that  there  is  rather  good  correlation 
between  the  static  and  dynamic  data. 

The  main  results  of  the  quasi-static  study  were: 

•  The  loading  isc therm  may  be  taken  as  unique,  independent  of 
intermediate  unload-reload  cycles. 

•  Unloading  and  reloading  occur  along  P-V  lines  with  moduli 
similar  to  those  of  the  solid  material. 

•  The  isotherms  show  an  initial  steep  loading  to  a  yield  point 
and  then  a  more  gradual  loading  up  toward  consolidation. 

Complete  consolidation  appears  to  be  approached  asymptotically. 

0  The  static  plain-strain  yield  point  is  between  50  and  100  percent 
of  the  Hugoniot  elastic  limit  and  the  correspondence  between  static 
and  dynamic  data  is  very  close  at  higher  stress  levels. 

2.  DYNAMIC  MEASUREMENTS  OF  SHOCK  ATTENUATION 

To  aid  in  development  and  verification  of  the  prediction  capability, 
several  shock-attenuation  experiments  were  performed  on  samples  of  porous 
copper,  iron,  and  tungsten  of  the  same  type  and  from  the  same  supplier  as 
specimens  used  by  Schmidt  and  Linde  (Ref.  5).  Specimens  were  cleaned  of 
pore  contaminants,  and  dynamic  experiments  were  performed  as  described  in 
Ref.  5.  Manganin  transducers  of  the  nominally  10-ohm  design  were  used  in 
all  experiments,  with  the  gage  element  located  nominally  0.25  mm  from  the 
foam  specimen  surface.  The  precise  gage  position  and  planarity  were 
neasured  optically  after  potting  of  the  gages  in  C-7  epoxy. 

Longitudinal  acoustic  velocities  were  measured  for  each  specimen 
prior  to  shock  loading.  These  measured  velocities  exhibited  sample-to- 
sample  variations  of  several  percent  and  were  generally  a  few  percent 
lower  than  the  velocities  measured  for  specimens  used  by  Schmidt  and  Linde 

(Ref.  5) 
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DATA  FROM  REF.  5 

O  INITIAL  DENSITY  6.31  TO  6.82g/cm3 

DATA  FROM  REF.  8 
•  INITIAL  DENSITY  =  6.052  g/cm  3 
A  INITIAL  DENSiTY  =  7.406g/cm3 
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HUGONIOT  FOR  SOLID  (REF.  5) 


r 


OTO  0.1 1  "  0J2  013  0JI4  oTi  OJ? 


SPECIFIC  VOLUME  - cm3/g 


FIG.  16  STATIC  AND  DYNAMIC  COMPRESSION  DATA  FOR  POROUS  COPPER 
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FIG.  17  STATIC  AND  DYNAMIC  COMPRESSION  FOR  POROUS  TUNGSTEN 


A  summary  of  the  shock  attenuation  experiments  is  given  in  Table  1. 

3.  STATIC  AND  DYNAMIC  COMPRESSION  OF  POL  YU  RE  THANE  FOAM 

The  backing  material  for  the  fly  r  plate  in  most  of  the  shock 
attenuation  tests  was  a  polyurethane  foam  with  a  Tensity  of  about 
0.185  g/cm3 .  The  Hugoniot  of  this  foam  was  estimated  through  a  series 
of  quasi-static  compression  tests  and  two  dynamic  experiments.  The 
quasi-static  compression  tests  were  performed  in  a  manner  similar  to 
that  used  for  the  porous  metals  except  that  the  engirdling  die  was 
much  lighter.  The  Hugoniot  measurements  (summarized  in  Table  II)  were 
made  with  a  copper-polyurethane  impact  and  a  polyurethane-polyurethane 
impact.  Two  gages  were  used  at  different  thicknesses  in  the  foam  on 
one  Hugoniot  shot  to  ascertain  that  a  steady  state  shock  front  had  been 
established.  The  results  of  both  static  and  dynamic  experiments  are 
shown  in  Fig.  18.  For  comparison,  the  solid  polyurethane  data  of 
Butcher  (Ref.  52)  is  also  exhibited. 


*  Actual  profiles  and  wave  arrival  times  obtained  are  given  in  Figs.  19 
through  24  of  Section  V. 
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Table  i 


CONDITIONS  Ci  ATTENUATION  EXPERIMENTS 


Porous 

Target 

Solid  Flyer 

Tilt 

Across 

Cage 

(Msec) 

Shot 

No. 

Matt  rial 

Thlcknear 

(cm) 

Density 
(g/cm3 ) 

Acoust  le¬ 
vel  ocl  ty 
cm/ sec 
( x  10“ ) 

Material 

Thickness 

(cm) 

Density 

(g/cm3) 

Vel ocl ty 
cm/sec 
(x  10‘) 

Iron 

0.297 

5.31 

3.10 

Iron 

0.0445 

7.81 

7.33 

<  0.114 

13,347 

0.158 

5.42 

3.31 

0.0465 

7.81 

7.42 

0.048 

13,403 

0.158 

5.29 

3.50 

0.165 

7.81 

7.38 

0.068 

13,418 

Copper 

0.314 

6.53 

3.4) 

Copper 

0.0325 

8.94 

7.35 

0.014 

13,3-19 

0.191 

6.43 

3.48 

0.0376 

8.94 

7.43 

0.054 

13,404 

0.159 

6.54 

3.32 

0.1651 

8  94 

7.43 

0. 140 

13,417 

Tungsten 

0.162 

14.56 

3.90 

Tungs  ten 

0.0505 

19.6 

7.38 

— 

13,402 

0.310 

14.26 

3.79 

0.114 

19.6 

7.35 

0.034 

13,473 

Notea 


1.  Basking  material  was  polyurethane  foam  >  0.183)  for  all  flyers  except  the  last. 

In  that  teat  the  backing  material  wja  Plexiglas  (Rohm  and  Haas  Co.)  with  an  initial  density  of  1 .  ’  89  g/cm3. 

2.  Hicorda  were  obtained  with  manganln  wire  gages  embedded  1  mil  into  C-7  epoxy.  Cages  were  immediately 
behind  the  target. 


Table  II 

HUGONIUT  DATA  FOR  POLYUR ETHiiNE  FOAM 


— 

Target 

Flyer 

Results 

Shot 

No. 

Thick. 

(mm) 

Density 
(g/cm3 ) 

Longlt. 

Acoustic 

Velocity 

(mm/»sec) 

Mat 1 1 . 

Thick. 

(mn.) 

Density 
(g/cm3 ) 

Impact 

Velocity 

rut 

Ac  rows 

Oage 

(k*sec) 

Stress 

(kbar) 

Particle 

Velocity 

(um/uaec) 

1 

* 

Denal ly 
(g/cm3) 

Shock 

Velocity 

(maiAisec) 

13405 

6.14 

0.186 

\ .  36 

** 

PF 

10.62 

0.187 

0.744 

0.031 

2.15 

0.366 

0.70 

0.50 

13420 

3.12 

C.  188 

1.55 

Cu 

6. 30 

8.94 

0.719 

- 

7.2 

0.702 

0.96 

0.87 

6.52 

7.2 

0.702 

*  Assuming  single  rhock. 
*♦  Polyurethene  foam. 
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SECTION  V 


COMPARISON  OF  THEORY  AND  EXPERIMENT 

1.  DISCUSSION  OF  WAVE- PROPAGATION  PHENOMENA 

Experimental  data  obtained  from  either  mechanical  impact  or  radiation 
deposition  experiments  generally  are  in  the  form  of  stress  or  particle 
velocity  histories  at  some  distance  from  the  impact  or  irradiated  surface. 

From  such  records,  we  can  estimate  or  assume  the  sequence  of  events  that 
resulted  in  the  record.  If  an  adequate  "equation  of  state"  and  wave- 
propagation  code  are  available  for  porous  materials,  such  histories  will 
be  predictable. 

For  impact  experiments,  it  is  desirable  to  use  a  solid  flyer  plate  of 
the  same  material  as  that  of  the  matrix  of  the  porous  target.  For  attenuation 
studies,  the  flyer  is  considerably  thinner  than  the  target.  According  to 
our  computations,  the  result  of  the  impact  is  often  a  series  of  impacts 
and  rebounds.  In  some  cases,  the  flyer  and  target  actually  separate;  but 
in  others,  the  stress  history  at  the  impact  interface  shows  only  a  series 
of  stress  oscillations,  each  oscillation  being  smaller  than  the  one 
preceding  it.  This  process  is  repeated  throughout  the  time  of  interest. 

To  follow  such  a  phenomenon  reliably  the  computer  code  should  allow  for 
spalling  and  recombination. 

In  the  case  where  the  porous  material  is  irradiated,  the  surface 
material  is  vaporized,  an  inner  section  is  melted,  and  the  deeper  material 
is  left  in  the  porous  state.  The  vapor  is  under  high  pressure  and  attempts 
to  expand  both  forward  and  backward.  The  forward  moving  vapor  is  termed 
the  blowoff.  The  molten  porous  material  is  like  a  mist  and  has  no  signif¬ 
icant  pressure  associated  with  it.  The  cooler  porous  material  is  under 
some  pressure  because  of  the  tendency  to  expand  under  heating.  The 
backward-moving  vapor  extends  rapidly  in  the  direction  of  the  molten 
material,  consolidating  it.  The  cooler  material  on  the  other  side  of  the 
molten  region  expands  much  more  slowly  toward  the  molten  material.  After 

some  time,  the  pressure  pulse  caused  by  vapor  expansion  meets  the 
expanding  cooler  material,  thus  completing  consolidation  of  the  molten 
material.  The  waves  that  travel  deep  into  the  porous  material  and 


reflect  from  the  rear  surface  are  the  initial  comp;  ession  stress  (caused 
by  radiative  heating),  then  an  expansion  wave  associated  with  the  motion 
of  the  cooler  material  toward  the  molten  material,  and  finall>  a  compressive 
wave  resulting  from  impact  of  the  vapor  (or  consolidated  molten  material) 
with  the  cooler  material.  This  third  wave  is  also  the  one  that  might  lie 
associated  with  the  blowoff  of  vaporized  material  from  the  front  surface. 
Because  of  the  dissipative  nature  of  the  porous  material,  the  third  wave 
will  usually  he  greatly  attenuated  before  it  reaches  the  rear  surface. 

2.  COMPUTED  STRESS  HISTORIES 

The  main  test  of  the  validity  of  both  the  computer  code  and  the 
equation  of  state  for  porous  materials  is  a  comparison  of  computed  stress 
histories  with  those  recorded  by  a  stress  gage.  The  experimental  situation 
of  most  interest  is  one  in  which  stresses  attenuated.  This  situation 
is  of  interest  both  because  of  the  intended  use  of  porous  materials  as 
a  stress-attenuating  medium  and  because  in  such  a  case  both  loading  and 
unloading  occur  so  that  virtually  all  aspects  of  the  theoretical  equation 
of  state  are  exercised. 

The  numerical  parameters  in  the  constitutive  relations  of  the  porous 
materials  were  selected  on  the  basis  of  experimental  Hugoniot  data  (Ref.  5) 
and  the  results  of  static  compression  tests.  These  basic  data  and  the 
selected  theoretical  isotherms  are  shown  in  Figs.  19  through  22.  The 
same  equation  of  state,  approximate  cell  sizes,  and  amounts  of  artificial 
viscosity  were  used  for  all  computations  for  a  material,  except  that  minor 
adjustments  (motion  of  the  curve  to  the  right  or  left)  were  made  in  the 
low-pressure  region  of  the  equation  of  s^te  to  account  for  the  slightly 
different  initial  density  of  each  sample. 

Computed  stress  histories  at  the  Manganin  gage  locations  and 
corresponding  experimental  gage  records  are  shown  in  Figs.  23  through  30. 
Stress  levels  are  those  attained  in  C-7  epoxy  about  0.25  mm  from  the 
specimen-epoxy  interface.  It  may  be  noted  that  the  most  detail  is  present 
in  the  records  in  which  a  great  deal  of  attenuation  has  occurred  and  there 
is  significant  separation  of  precursor  and  main  wave.  These  records  are 
most  instructive  for  verifying  the  constitutive  relations  of  the  porous 
material . 
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DATA  FROM  REF.  5 
(INITIAL  DENSITY  =  5.5g/cm3) 

0  AS  -RECEIVED  MATERIAL 
A  CLEANEO  MATERIAL 


DATA  FROM  REF.  II 
•  INITIAL  DENSITY  *  4.5g/cm3 
A  INITIAL  OENSITY=  5.8g/em3 
■  INITIAL  DENSITY  =  7.0  g/cm5 


HUGONIOT  FOR  SOLID  (REF  5) 


P-V  RELATION  FOR  COMPUTATIONS 


STATIC  OATA 


**  o 

^  A  <5*. 


•  YIELD 


SPECIFIC  VOLUME - cmVg 


FIG.  19  COMPARISON  OF  THEORETICAL  ISOTHERM  WITH  DATA  FOR  POROUS  IRON 
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DATA  FROM  REF.  5 

O  INITIAL  DENSITY  6.31  TO  6  82  g/cm5 

DATA  FROM  REF.  8 
«  INITIAL  DENSITY  =  6.052g/cmS 
A  INITIAL  DENSITY1  7  406g/cm3 


HUGONIOT  FOR  SOLID  (REF. 5) 


\  /  P-V  RELATION  FOR  COMPUTATIONS 


,  STATIC  DATA 


YIELD  A 


YIELD 


SPECIFIC  VOLUME  -  cm5/g 


FIG.  20  COMPARISON  OF  THEORETICAL  ISOTHERM  WITH  DATA  FOR  POROUS  COPPER 


FIG.  21  COMPARISON  OF  THEORETICAL  ISOTHERM  WITH  DATA 
FOR  POROUS  TUNGSTEN 
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FIG.  22  COMPARISON  OF  THEORETICAL  ISOTHERM  WITH  DATA  FOR 
POLYURETHANE  FOAM 
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FIG.  25  IRON/POROUS  IRON  IMPACT.  SHOT  13418 
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TIME  - ^**c 

BA-65S6-57 


FIG.  29  TUNGSTEN/POROUS  TUNGSTEN  IMPACT,  SHOT  13402 
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The  comparisons  for  porous  iron,  Figs.  23  through  25  show  fairly 
good  correlation  of  wave  front  shapes;  the  precursors  have  the  right 
amplitudes  and  even  tht  precursor  reverberations  (small  jumps  in  stress 
between  the  precursor  and  main  wave)  are  apparent  in  Fig.  23.  The 
computed  magnitude  of  the  main  wave  is  an  overestimate  in  Figs.  23  and 
25,  but  is  correct  in  Fig.  24.  The  arrival  times  coincide  satisfactorily 
in  Figs.  23  and  25,  but  not  in  Fig.  2a.  The  variability  of  the  correla¬ 
tion  is  inveresting  because  the  same  equation  of  state  was  used  for  all 
of  the  computations.  The  most  probable  causes  for  the  unevenness  in  the 
correlation  are: 


•  Dominance  of  different  features  of  the  equation  of  state  in  the 
three  cases;  e.g.,  the  high-pressure  portion  of  the  Hugoniot  and 
the  unloading  modulus  are  most  important  for  a  thin  target 
under  high-velocity  impact,  while  the  low-pressure  portion 

of  the  Hugoniot  is  most  important  if  the  stress  at  the  gage  is 
low.  However,  this  reason  is  probably  not  important  for  the  iron 
impacts;  the  test  conditions  of  Figs. 23  and  25  bracket  those  for 
Fig.  24. 

•  Material  variability  among  the  three  samples  tested.  For  example, 
even  though  samples  were  "cleaned"  (Ref.  5),  the  measured 
acoustic  velocities  (  see  Table  I)  are  somewhat  more  variable 
than  those  obtained  earlier  on  similar  material  (see  Ref.  5), 

and  even  the  earlier  specimenss  showed  appreciable  sample- 
to-sample  variation. 

The  comparison  of  stress  histories  for  the  porous  copper  impacts 
are  shown  in  Figs.  26  through  28.  The  computed  arrivals  are  a  little 
late  in  Figs.  26  and  27,  but  the  peak  stresses  compare  satisfactorily 
with  the  experimental  records.  The  arrival  time  is  satisfactory  in 
Fig.  28  but  the  ccmputed  peak  stress  is  much  higher  than  the  measured 
peak  stress.  The  marked  peaks  in  both  experimental  and  computed  records 
correspond  to  actual  reverberations  through  the  flyer  plate.  The  computed 
results  for  this  case  are  strongly  dependent  on  the  thickness  of  the  flyer 
plate  and  on  the  unloading  wave  velocity.  The  stress  at  the  point  of 
impact  was  70  kbar.  This  high  stress  persists  almost  to  the  gage.  Then 
the  first  unloading  wave  arrives  and  reduces  the  stress  at  the  gage. 

With  a  5  percent  increase  in  flyer  plate  thickness  or  5  percent  reduction 
in  unloading  velocity,  the  peak  stress  at  the  gage  would  have  been  70  t  bar. 
In  addition  to  the  possible  sources  of  discrepancy  discussed  above  ror 
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iron,  possible  causes  for  the  disagreement  between  computed  and  measured 
stress  histories  of  Fig.  28  are  as  follows: 

•  The  high  tilt  that  occurred  in  the  experiment  (see  Table  I) 
would  lead  to  an  apparent  damping  of  measured  peaks. 

•  An  error  may  have  occurred  in  measuring  the  flyer  thickness. 

•  The  solid  equation  of  state  may  be  too  "soft,''  hence  velocities 
are  too  low.  A  higher  velocity  would  not  significantly  alter 
the  results  of  Figs.  26  and  27  but  could  drastically  alter  the 
correlation  in  Fig.  28. 

Clearly,  a  repeat  of  Shot  13,417  should  be  performed. 

Two  tungsten  records  are  shown  in  Figs.  29  and  30.  The  correlation 
between  measured  record  and  computed  history  for  both  stress  amplitudes 
and  arrival  times  appears  to  be  acceptable. 

The  computed  results  are  primarily  a  function  of  the  experimental 
conditions  and  the  theoretical  Hugoniot  but  are  also  modified  to  some 
extent  by  the  cell  size,  amount  of  artificial  viscosity,  the  time  and 
extent  of  rezoning,  and  other  discretionary  factors.  Therefore,  for  the 
convenience  of  those  wishing  to  reproduce  these  records,  the  complete 
data  input  for  each  of  the  eight  computations  is  given  in  the  Appendix. 

The  results  depicted  in  the  previous  eight  figures  and  the  calcula¬ 
tions  that  lead  to  these  as  a  final  result  have  brought  the  authors  to 
the  following  two  main  conclusions: 

1.  The  constitutive  relations  (or  equation  of  state)  and  the 
computational  scheme  developed  for  porous  materials  are 
adequate  for  predicting  stress  wave  magnitudes  and  profiles 
in  a  porous  material  as  a  function  of  time.  The  precursor, 
small  reflected  waves,  and  other  details  of  wave  propagation 
in  a  porous  material  are  adequately  represented;  hence,  the 
model  may  be  used  to  study  stress  wave  interactions,  including 
details  of  a  compressive  wave  front  and  subsequent  rarefactions. 
The  major  uncertainty  is  the  numerical  values  of  parameters  for 
a  given  material  having  specimen-to-specimen  variations. 
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2.  The  uncertainty  in  the  results  obtained  from  discretionary 
factors  in  the  present  computation  scheme  can  be  as  great  as 
the  difference  between  the  experimental  and  computed  records. 

This  uncertainty  can  be  examined  by  doubling  or  halving  the 
cell  size  and  artificial  viscosity  values  that  were  used.  Such 
changes  will  modify  peak  stresses  by  5  to  10  percent  and  alter 
arrival  times  by  smaller  percentages.  Although  neither  cell 
size  nor  artificial  viscosity  has  an  intended  physical  signifi¬ 
cance,  its  effect  on  the  computed  results  is  as  real  as  are 
changes  in  the  theoretical  Hugoniot.  The  ideal  value  of  cell 
size  is  the  minimum  that  can  be  afforded  (computation  times 
increase  as  the  inverse  square  of  cell  size).  Appropriate 
values  of  artificial  viscosity  are  suggested  in  the  next  section. 

3.  CONSTRUCTION  OF  INPUT  DATA 

The  details  of  writing  the  input  for  the  computer  program  are  given 
in  the  manual  for  the  SRI  PUFF  1  code  (Ref.  50).  in  this  section  some 
guidance  for  developing  input  quantities  is  presented,  based  on  our 
experience  with  the  program. 

The  first  step  is  to  select  a  realistic  Hugoniot  (actually  a 
'dynamic  isotherm '  is  used  in  the  code)  for  the  material .  As  evidenced 
in  Figs.  15  through  18  both  static  and  dynamic  data  are  useful  in  con¬ 
structing  tnls  theoretical  Hugoniot.  The  Hugoniot  need  not  be  smooth: 
even  severe  Joints  lead  to  only  minor  oscillations  in  the  stress  records 
at  early  times  in  the  computation.  The  EMELT  variables,  which  define 
the  variation  of  strength  with  internal  energy,  have  been  selected  by  using 
published  melting  point  and  enthalpy  data  plus  information  on  the  variation 
of  yield  or  ultimate  strength  with  temperature. 

The  initial  bulk  and  shear  moduli  a*  e  selected  based  on  the  known 
or  expected  precursor  velocity.  The  slope  of  the  initial  section  of 
Hugoniot  is  also  constructed  to  correspond  with  the  precursor  velocity. 

The  strain  hardening  parameter,  YADDP,  is  varied  through  the  porous 
regions  to  provide  for  an  increase  of  yield  strength  so  that  the  yield 
is  that  of  the  solid  when  the  material  is  compacted  to  solid. 


The  layout  of  cells  for  an  impact  problem  is  made  with  small  cells 
near  the  point  of  impact  and  large  cells  farther  away.  It  was  found 
that  to  achieve  fidelity  to  detail  for  the  impact  problems,  at  least  10 
cells  should  be  used  for  a  thin  solid  flyer.  If  the  wave  front  is  of 
paramount  importance,  then  more  cells  should  be  used  throughout.  For 
impact  conditions  similar  to  tnose  studied  here,  the  cells  in  the  porous 
target  at  the  impact  point  should  be  about  half  the  size  of  the  cells 
in  the  flyer.  This  is  in  accord  with  the  common  practice  of  matching 
travel  time  through  cells  at  interfaces. 

These  small  cells  are  necessary  because  the  peak  stress  at  the 
interface  usually  is  not  reached  until  the  first  porous  cell  is  compacted 
to  solid.  Because  of  the  motions  required  to  produce  compaction,  some 
time  is  needed  to  consolidate  the  first  cell  and  ,.s  time  appears  as 
the  initial  rise  of  the  stress  wave.  If  the  initial  rise  is  comparable 
to  the  duration  of  the  stress  wave  (twice  the  travel  time  through  the 
flyer),  then  the  interface  stress  will  not  resemble  the  "true"  stress 
wave. 

Even  if  the  impact  is  instantaneous,  the  resulting  stress  wave  will 
broaden  after  it  has  propagated  some  distance  into  the  porous  material. 
Hence,  at  later  times  and  greater  depths  into  the  material  the  cells  need 
not  be  so  small.  Our  practice  has  been  to  initiate  the  problem  with  small 
cells  at  the  impact  point  in  the  target,  varying  up  to  cells  3  to  10  times 
as  large  at  the  rear  of  the  target.  Then  at  the  conclusion  of  the  main 
impact  with  the  flyer  (one  reverberation  of  stress  through  the  flyer),  the 
REZONE  subroutine  is  called  and  the  small  cells  near  the  impact  point  are 
increased  several  fold  so  that  all  cells  in  the  target  are  about  the  same 
size. 

For  radiation  deposition  problems  the  cells  should  be  small  at  the 
irradiated  surface  and  may  increase  in  size  in  the  cooler  region.  At 
least  4  or  5  cells  should  be  in  the  region  thee  is  vaporized.  If 
possible  economically,  small  cells  should  be  used  also  throughout  the 
molten  zone  to  minimize  oscillations  in  the  calculated  stress.  The  cells 
may  be  rezoned  following  completion  of  deposition,  but  only  a  small 
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increase  (e.g.,  doubling)  in  cell  sizes  should  be  permitted  at  that  time. 
Later,  after  the  molten  material  has  been  compacted  and  the  vapor  is  well 
expanded,  the  cells  may  be  rezoned  to  several  times  their  original  size. 

For  both  impact  and  radiation  problems  the  artificial  viscosity  may 

* 

be  set  at  COSQ  -  4.0  and  Cl  =0.1  for  the  first  try.  Then  a  computation 
should  be  run  and  the  results  examined.  The  results  may  indicate  that 
changes  in  the  input  data  are  desirable.  The  following  are  guidelines 
in  d-?ciding  which  data  to  change: 

•  If  your  Hugoniot  is  not  well  established,  then  note  that  increases 
in  the  concavity  of  the  Hugoniot  will  increase  stress  attenuation 
and  will  usually  cause  the  main  wave  to  arrive  a  little  later 
(because  lower  stresses  will  be  associated  with  lower  wave  velocities). 

•  A  decrease  in  the  cell  dimensions  makes  a  sharper  delineation  of 
the  wave  front,  showing  yield  point  and  reflected  waves  more 
clearly.  Smaller  cells  also  reduce  the  stress  attenuation, 
steepen  wave  fronts,  and  may  augment  oscillations  in  the  wave  front. 

•  A  decrease  in  the  artificial  viscosity  may  induce  instability 

but  will  usually  only  allow  an  increase  in  oscillations.  However, 
the  decrease  may  also  decrease  oscillations  if  there  had  been  too 
much  viscosity  before.  The  decrease  in  viscosity  will  also  de¬ 
crease  the  attenuation,  steepen  wave  fronts,  sharpen  the  definition 
of  features  in  the  wave  front,  and  retard  arrival  times  (viscosity 
slightly  augments  wave  velocity). 


*  For  steady-state  wave  propagation,  the  physically  correct  amount  of 
viscosity  is  that  which  causes  the  total  stress  state  to  follow  a 
straight  line  from  yield  to  the  peak  stress  on  a  stress-volume  plot 
(Rayleigh  line) .  The  correctness  can  only  be  judged  after  a 
computation  by  plotting  successive  values  of  R  (total  stress)  versus 
1/p  (p  =  density)  for  a  few  cells.  (Cells  near  an  impact  interface 
should  not  be  considered  in  the  comparison  because  their  stress- 
volume  paths  will  not  follow  the  Rayleigh  line.) 


4.  AREAS  REQUIRING  FURTHER  EFFORT 

During  the  construction  of  the  equation  of  state  and  the  computer 
program  we  discovered  some  problem  areas  outside  tne  scope  of  the  present 
work.  The  most  notable  problem  was  the  lack  of  a  complete  spall  criterion. 
Only  a  rudimentary  one  has  been  inserted  during  the  project  but  this  is 
not  satisfactory  for  several  reasons.  The  criterion  is  independent  of 
time  although  it  is  known  that  spall  strength  is  a  function  of  the 
tensile  stress  duration.  The  rudimentary  criterion  is  also  independent 
of  internal  energy  so  that  a  vapor  has  the  same  spall  strength  as  the 
solid.  Some  of  these  more  glaring  deficiencies  could  be  easily  remedied 
in  the  program  but  it  was  felt  that  such  changes  should  await  the  results 
to  be  obtained  under  a  current  project  at  SRI  on  fracturing. 

In  code  computations  of  radiation  deposition  in  solid  materials, 
there  are  often  violent  nonphysical  oscillations  of  pressure  in  the 
vaporized  region.  When  the  material  is  initially  porous,  these 
oscillations  also  occur,  possibly  with  even  more  violence.  For  porous 
materials  one  cause  of  the  oscillations  may  be  the  discontinuities  in  the 
theoretical  formulation  of  the  equation  of  state  surface.  Such 
discontinuity  occurs  between  the  zero  pressure  plane  of  the  molten  mist 
and  the  surface  for  the  expanding  vapor.  When  a  cell  passes  through 
this  discontinuity  (when  the  material  is  consolidated  from  a  distended 
mist  to  a  poreless  liquid  or  gas),  large  stress  oscillations  are  induced 
in  the  consolidated  material.  Such  oscillations  are  nonphysical  and 
corrections  should  be  made  in  the  program  to  remove  them.  The  oscillations 
are  of  such  a  magnitude  that  they  make  nonsense  of  spalling  calculations 
in  the  molten  and  vaporized  material.  Because  the  motion  of  the  molten 
and  vaporized  material  has  an  important  effect  on  the  blow-off  momentum, 
and  therefore  on  the  momentum  transmitted  through  the  porous  material 
to  the  backing  material,  stress  oscillations  are  of  great  concern.  The 
consolidation  that  induces  the  oscillations  is  similar  to  the  situation 
occurring  during  an  impact  and  could  be  handled  with  similar  steps. 

•  Use  a  special  iterative  interface  computation  at  the  boundary 
between  vapor  and  molten  material. 
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Employ  very  short  time  steps  for  20  cycles  just  following 
consolidation. 


•  Increase  the  viscosity  of  the  interface  between  vapor  and 
molten  material. 

•  Use  smaller  cells  in  the  consolidating  region. 

Another  change  that  would  reduce  the  severity  of  the  impact  and  thus 
reduce  the  magnitude  of  the  oscillations  is  to  provide  for  the  partial 
pressure  of  the  molten  material,  i.e.,  make  an  adjustment  in  the 
expansion  equation  of  state.  Some  combination  of  these  approaches  should 
be  used  to  improve  the  calculation  of  the  motion  of  the  vapor. 

In  experimental  studies  of  porous  materials,  a  large  number  of 
repetitive  tests  should  be  included.  There  appears  to  be  much  more 
scatter  in  experimental  data  from  distended  material  than  from  the 
corresponding  solids.  This  scatter  is  to  be  expected  because  tne 
constitutive  relations  may  depend  on  the  pore  size,  previous  stress 
history  of  the  particles,  interparticle  bonding,  inclusions,  and  other 
factors  that  may  accidentally  vary  from  sample  to  sample.  Therefore, 
instead  of  expecting  single  values  for  data  points,  we  should  obtain 
mean  values  and  standard  deviations  from  the  values.  Then  our  mathe¬ 
matical  model  can  be  constructed  to  represent  the  average  behavior  of  the 
material.  As  a  basis  for  such  an  average  model,  it  is  necessary  to 
repeat  experiments  under  apparently  identical  conditions.  For  basic 
studies  on  porous  materials  great  care  should  be  taken  to  acquire 
samples  with  uniform  properties.  For  Hugoniot  and  attenuation 
experiments,  a  minimum  of  three  identical  shots  should  be  conducted. 

We  should  start  by  repeating  the  shock  attenuation  experiments  reported 
above  (particularly  Shot  13,417). 

The  constitutive  relations  developed  on  the  project  have  been 
shown  to  represent  the  response  of  a  room-temperature  porous  material 
undergoing  a  low-velocity  impact.  Such  an  impact  exercises  the 
constitutive  relations  only  in  the  vicinity  of  the  Hugoniot.  To  more 
completely  verify  these  relations,  it  is  necessary  to  study  high- 
pressure  (megabar)  impacts,  externally  heated  specimens,  and  the  case 
of  radiation  deposition.  Sparse  experimental  data  of  both  types  are 


available.  More  data  should  be  obtained  and  compared  with  shock^-wave 
computations,  using  the  present  constitutive  relations  to  further  verify 
or  modify  the  relations. 

A  great  deal  of  intuition  and  very  little  experimental  evidence 
has  been  employed  to  define  the  detailed  behavior  of  the  porous  material 
during  loading  and  unloading,  heating,  cooling,  yielding,  and  spalling. 

As  a  consequence,  there  are  many  small  uncertainties  in  the  formulations. 
Some  questions  that  have  not  yet  been  treated  in  detail  but  could  readily 
be  further  investigated  are: 

•  How  long  does  it  take  for  the  porous  material  to  come  to 
thermal  equil  ibrium?  To  what  degree  of  accuracy  can  the 
"consolidated”  material  be  represented  by  the  usual  equation 
of  state  for  .•*.  solid?  Are  some  characteristics  of  the 
uncompacted  material  retained? 

•  What  are  the  microscopic  details  of  the  collapse  of  the 
pores.  How  great  are  the  thermal  gradients  within  the 
yielding  particles?  How  can  we  theoretically  derive 
appropriate  values  for  the  "effective"  Gruneisen  ratio? 

•  What  is  the  effect  of  solid,  liquid,  or  gaseous  inclusions 

in  the  pores?  Such  inclusions  may  have  importance  in  allowing 
foams  to  operate  effectively  as  countermeasure  materials  in 
"shine- through"  situations  where  premature  collapse  of  the 
foams  would  render  them  ineffective. 

Some  of  these  problems,  such  as  stress  oscillations  in  the  vapor, 
temperature  effects  and  off-Hugoniot  behavior,  and  the  spalling  criterion, 
should  be  handled  immediately;  the  other  problems  may  involve  refine¬ 
ments  that  are  not  economically  justifiable  for  present  '.pplications . 
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APPENDIX 


INPUT  DATA  FOR  COMPUTED  STRESS  RECORDS 

The  SRI  PUFF  1  computer  code  was  employed  to  calculate  the  stress 
records  shown  jn  Figs.  23  through  30.  The  Input  data  on  which  the 
calculations  were  based  are  given  in  the  following  pages. 
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STRESS  HISTORY  AT  INTERFACES  BETWEEN  MATERIALS  1  AND  2  IS  FOP  COMPARISON  WITH 
CAGE  RECORD  OF  OAVE  SCHMIDT  FROM  SHOT  13147 

COMPUTATIONS  WERE  MAOE  WITH  SRI  PUFF  1  ON  THE  COC  3200  AT  SRI.  THF  VFRSION 
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NPOR  =  0 

NYO  = 

n 

NCON  =  0 

EQSTC  * 

7.816EM0 

EOSTO  = 

1 .956EU0 

EQSTE  = 

3.000E+1 0 

EOSTO 

= 

2.0U0F+00 

EQSTH  * 

2. 500E-01 

EOSTS  = 

2.2  1 4E  -*■  1 1 

COSQ  * 

2.000E+0C 

Cl  - 

0.050E  00 

C2  = 

0 

TENS! 11*  - 

3. OOOE+09 

TENS! 2  5* 

0. 

TENSI3 )= 

0. 

NZONES*  It 

50 

CELLS  IN 

1 . 5006*00 

CM,  DELX= 

1 . 200E-02 

*FND* 

80 


raK?’' 

'  h 


NDATE*  10/22/68  I0ENT-1O44  IRON-IRON  FOAM  IMPACT  13403 

STRESS  HISTORY  AT  INTERFACES  8ETWEFN  MATERIALS  1  AND  2  IS  FOR  COMPARISON  W1TH 
GAGE  RECORD  OF  DAVf  SCHMIDT  FROM  SHOT  13403 

COMPUTATIONS  WERE  MADE  WITH  SRI  PUFF  1  ON  THE  COC  ?200  AT  SRI.  THF  VERSION 
OF  SRI  PUFF  1  IS  EQUIVALENT  TO  THAT  SUPPLIED  TO  AFWL  WITH  FINAL  RFPOPT  ON 


DISTENDED 

MATERIAL  1 

MODEL  DEVELOPMENT. 

* 

*  *  IRON- 

-IRON  FOAM 

IMPACT  * 

*  4 

1  NTEPT  * 

6 

NJEDIT  * 

6 

NREZON  * 

1 

NSEPRAT  = 

0 

2  TEDI TS* 

3.000E-07 

5. OOOE-O  7 

1 .OOOE-06 

1. 500E-0& 

2.000F-06 

2. 500E-06 

3.000F-06 

3  JED  ITS' 

28  33 

38  43 

48  56 

4  NTR  * 

1 

5  JREZON* 

48 

6  NEDTM  - 

10000 

NED  IT  * 

10003 

NPERN  * 

3 

7  STOPS 

JCYCS  * 

1000 

CKS  * 

3.000F  on 

TS  = 

3.0006-06 

8  NMTRLS* 

4 

MAl'FL  * 

2 

UZERO  * 

7.420E  04 

POLYURETHANE  FM 

RHOS  = 

1.19 

NEQST  -  1 

NPOR  *  4 

NYO  =  1 

NCON  *  0 

EQSTC  - 

7.816E*10 

EQSTO  = 

1.956E*11 

EOSTE  * 

3.0006*10 

EOSTC.  = 

?.000F*00 

EOSTH  * 

0.250E  *00 

EQSTS  = 

2.2146*11 

YOS  = 

1.0006*00 

MU  = 

1.0006*09 

YADD  * 

0. 

COSQ  = 

2.0006*00 

Cl  * 

0.050E  00 

C2  * 

0 

TENS!  1)*  - 

-l.OOOEMO 

TENS! 2 ) =  - 

- 1 . OOOE* 10 

TENSI31*  - 

-l.OOOE+lO 

RHOP  = 

0.187 

0.202 

0.3  33 

0.625 

1.205 

1.250 

COSQ* 

2.0 

2.0 

2.0 

2.0 

2.0 

Cl* 

0.05 

0.05 

1.05 

0.05 

0.05 

1  P2  * 

2.8006*07 

DELP  * 

6. OOOE *06 

YADD  - 

0. 

2  P2  * 

3. 350E*07 

OELP  *  - 

- 1 ,400£*06 

YADD  - 

0. 

3  P2  * 

3.450E+07 

DELP  = 

2.500E+05 

ADD  * 

0. 

4  P2  * 

8.000E+08 

OCLP  *  - 

- 1.9006*08 

YADD  = 

0. 

EMELT  * 

3.0006*09 

9.0006*08 

0. 

0.5 

0. 

AK  - 

1.000E*09 

MUP  * 

1 .000E+09 

YO  * 

0. 

NZONES*  1. 

10 

CELLS  IN 

6.0006-01 

CM, 

DELE  IN* 

4.00OF-02 

IRON 

RHO  - 

7. B10E+00 

NEQST  *  l 

NPOR  *  0 

NYO  =  1 

NCON  *  0 

EQSTC  * 

1 . 568E* 12 

6QSTD  = 

5.644E+12 

EQS7E  * 

8.4006*10 

FQSTG  * 

1.600 

EQSTH  * 

0. 250E*00 

EOSTS  - 

1.2166*13 

YO  * 

5.741E>09 

MU  * 

6 . OOOE ♦ l 1 

YADD  = 

0. 

COSQ  * 

2.0006*00 

Cl  * 

0.050E  00 

C2  * 

0 

TENS! II*  - 

■l.OOOE* 1 1 

TENSI2I* 

0. 

TENS!  3  )*  -l.OOOEMO 

NZONES*  1. 

10 

CELLS  IN 

4. 650E-02 

CM 

IRON  FOAM 

RHOS  - 

7.81 

NEQST  -  l 

NPOR  *  4 

NYO  -  1 

NCON  «  0 

EQSTC  * 

1.5686*12 

EOSTD  * 

5.6446*12 

EQSTE  - 

8.400F*10 

FOSTG  * 

1.600 

EQSTH  * 

0. 2506*00 

EQSTS  * 

1.2 16E*13 

YO  * 

5«.  700E*09 

MU  - 

6. OOOE  *11 

YADD  = 

0. 

COSQ  * 

4.000E*00 

Cl  - 

0.100E  00 

C  2  * 

0 

TENS  III*  - 

1.0006*1 1 

TENS! 2 ) *  - 

■l .  OOOE*  10 

TENSC3I*  - 

-3.0006*09 

RHOP  * 

5.42 

5.442 

7.23 

7.73 

7.95 

8.0 

COSQ* 

4.0 

4.0 

4.0 

4.0 

4.0 

Cl* 

0.1 

0.1 

0.1 

0.1 

0.1 

l  P2  * 

1 .  800E*09 

DELP  * 

0. 

YADD  * 

0. 

2  P2  * 

8. OOOE+09 

DELP  -  - 

- 1 .OOOE  09 

YADD  * 

0. 

3  P2  * 

1.2006  10 

OFLP  =  - 

■  1 . 000E*09 

YADD  * 

0 

4  P2  * 

3.5006*10 

DELP  =  - 

•2. OOOE*09 

YADD  - 

3. ROOF  *09 

EMELT  = 

4.7506*09 

2 . OOOE  *09 

0. 

0.5 

0. 

AK  * 

6.0006*11 

MUP  * 

2 .OOOE  *  1 1 

YO  * 

0. 

NZONES*  1, 

30 

CELLS  IN 

1 . 580E-01 

CM,  DELX* 

2. 3006-03 

C 7  (MN  GAGE) 

RHO  * 

1.1906*00 

NEQST  *  l 

NPOR  *  0 

NYO  *  0 

NCON  «  0 

EQSTC  = 

7.R16E+10 

EQSTO  - 

1 • 956E  *10 

EQSTE  = 

3.0006*10 

FOSTG  * 

2.000F+00 

EQSTH  * 

2. 500E-0 1 

EOSTS  * 

?.2l4E*ll 

COSQ  = 

2.0006*00 

Cl  * 

0.050E  00 

C2  - 

0 

TENS  III*  - 

3.000F*09 

TENS!’)* 

0. 

TENS! 31* 

0. 

NZONES-  1. 

30 

CELLS  IN 

7.5006-01 

CM,  DELX* 

1 . 200E-02 

♦  END* 


81 
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NDfcTE*  10/22/68  IDENT*1045  IRON-IRON  FOAM  IMPACT  13*18 

STRESS  HISTORY  AT  INTERFACFS  BETWEEN  MATERIALS  1  ANO  2  IS  FOR  COMPARISON  WITH 
GAGE  RECORD  OF  DAVE  SCHMIDT  FROM  SHOT  13418 

COMPUTATIONS  WERE  MADE  WITH  SRI  PUFF  1  ON  THE  CDC  3200  AT  SRI.  THE  VERSION 
OF  SRI  PUFF  1  IS  EQUIVALENT  TO  THAT  SUPPLIED  TO  AFWL  WITH  FINAL  RFPORT  ON 


DISTENDED 

MATERIAL  1 

MODEL  DEVELOPMENT. 

* 

*  *  IRON- 

-IRON  FOAM 

IMPACT  ♦ 

*  * 

1  NTEDT  * 

4 

NJEDIT  * 

6 

NREZON  * 

1 

NSFPRAT 

S 

0 

2  TEDITS* 

1.000E-06 

1.500E-06 

2.000E-06 

2 • 500E-06 

3  JEDITS* 

58  63 

68  73 

78  86 

4  NTR  * 

1 

5  JREZCN- 

78 

6  NEDTM  x 

10000 

NEDIT  x 

10000 

NPERN  * 

3 

7  STOPS 

JCVCS  x 

1000 

CKS  * 

3.000F  00 

TS  - 

2 . OOOE-06 

8  NMTRLS* 

4 

MATFL  » 

2 

UZERO  * 

7.380E  04 

POLYURETHANE  FM 

RHOS  * 

1.19 

NEQST  =  l 

NPOR  =  4 

NYO  * 

1 

NCON  -  0 

EQSTC  « 

7.8 16E+10 

EQSTO  * 

1 .956E  *11 

EQSTE  * 

3. OOOE* 1 0 

EOS TO  * 

2. OOOF  *00 

EQSTH  * 

0. 250E+00 

EQSTS  * 

2.214E*11 

YOS  - 

1 . OOOE+OO 

MU  = 

1 . OOOE  *09 

YADD  * 

0. 

COSQ  - 

2.000E+00 

Cl  * 

0.050E  00 

C  2  * 

0 

TENS  I  l )»  - 

1.000E*10 

TENSI2I*  - 

■1.000E*10 

TENS (31*  ■ 

-1.000F+10 

RHOP  x 

0.187 

0.202 

0.333 

0.625 

1.205 

1.250 

COSQ* 

2.0 

2.0 

2.0 

2.0 

2.0 

Cl- 

0.0$ 

0.05 

0.05 

0.05 

0.05 

1  P2  x 

2.800E*07 

OELP  - 

6. OOOE  *06 

YADD  * 

0. 

2  P2  * 

3.350E*07 

DELP  x  - 

-1.400E*06 

YADD  * 

0. 

3  PI  * 

3.450E*07 

OELP  * 

2.500E  *05 

YADD  - 

0. 

4  P2  x 

8. OOOE *08 

DELP  -  - 

■1.900E*08 

YADD  * 

0. 

EMELT  « 

3.000E*09 

9.000E+08 

0. 

0.5 

0. 

AK  - 

1.000E*09 

MUP  » 

1.000E+09 

YO  - 

0. 

NZONES*  1, 

10 

CELLS  IN 

6.000E-01 

CM, 

DELFIN* 

4. OOOE-02 

IRON 

RHO  * 

7.810E+00 

NFQST  •  l 

NPOR  *  0 

NYO  * 

1 

NCON  =  0 

EQSTC  « 

1. 568E* 12 

EQSTO  * 

5 .644E  *  12 

EQSTE  * 

B.400F+10 

EQSTO  = 

1.600 

EQSTH  « 

0. 250E  *00 

EQSTS  * 

1.216E*13 

YO  * 

5.741E*09 

MU  * 

6. OOOE* 1 1 

YADD  » 

0. 

COSQ  x 

2.000E*00 

Cl  - 

0.050E  00 

C2  - 

0 

TENS 1 1 )x  - 1  .OOOE* 1 1 

TENSI2I- 

0. 

TENS! 3 )*  -1.000E+10 

NZONES*  1. 

40 

CFLLS  IN 

1.650E-01 

CM 

IRON  FOAM 

RHOS  * 

7.81 

NEQST  *  1 

NPOR  *  4 

NYO  * 

1 

NCON  * 

EQSTC  * 

I.S68E*12 

EQSTO  * 

5.644E*12 

EQSTE  * 

8.400E*1 0 

EOSTO  = 

1.600 

EQSTH  « 

0. 250E*00 

EQSTS  * 

1.216E*13 

YO  * 

5.700E+09 

MU  * 

6. OOOE  *11 

YADD  x 

0. 

COSQ  * 

4. 000E*00 

Cl  * 

O.IOOE  00 

C2  * 

0 

TENS! 1 )*  - 

1.000E*ll 

TENSI2I*  - 

•1.000E*10 

TENS! 3 )*  - 

-3.000E+09 

RHOP  x 

5.29 

5.312 

7.23 

7.73 

7.95 

8.0 

COSQ= 

4.0 

4.0 

4.0 

4.0 

4.0 

Clx 

0.1 

0.1 

0.1 

0.1 

0.1 

1  P2  x 

1.800E+09 

OELP  * 

0. 

YADD  - 

0. 

2  P2  x 

8.000E*09 

DELP  *  - 

■l.OOOE  09 

YADD  * 

0. 

3  PI  * 

1.200E  10 

DELP  «  - 

■l . 000E*09 

YADD  » 

0 

4  P2  * 

3. 500E*10 

DELP  *  - 

•2.000E+09 

YADO  * 

3. 800F+09 

EMELT  * 

4.750E+09 

2.000E+09 

0. 

0.5 

0. 

AK  x 

6.000E+11 

MUP  = 

2. OOOE* 1 1 

YO  - 

0. 

NZONES*  1, 

30 

CELLS  IN 

1  • 580E-01 

CM,  DELX* 

2.300E-03 

C7  I MN  GAGE) 

RHO  = 

1.1 90E  *00 

NEQST  =  1 

NPOR  *  0 

NYO  = 

0 

NCON  =  0 

EQSTC  * 

7.816E*10 

EOSTO  * 

1.956E+10 

EQSTE  * 

3. OOOE  *10 

EOSTO  = 

2.000E*00 

EQSTH  * 

2. 500E-01 

EQSTS  * 

2.2l4E*ll 

COSQ  x 

2.000E*00 

Cl  = 

0.050E  00 

C2  = 

0 

TENS  1 1 ) x  - 

3. 000E*09 

T  ENS (21  = 

0. 

TENS  1 3 )  = 

0. 

NZONES*  1« 

30 

CELLS  IN 

7.  500E-01 

CM,  DELX= 

1.200E-0? 

•  FNO* 


82 


NDATE*  10/17/68 


I0ENT-10A0  COPPER-COPPER  FOAM  IMPACT 


11141 


STRESS  HISTORY  AT  INTERFACES  BETWEEN  MATERIALS  1  AND  2  IS  FOR  COMPARISON  WITH 
GAGE  RECORD  OF  DAVE  SCHMtOT  FROM  SHOT  133A9 

COMPUTATIONS  WERE  MADE  WITH  SR!  PUFF  1  ON  THE  CDC  3200  AT  SRI,  ThF  VFRSION 
CF  SR!  PUFF  1  IS  EQUIVALENT  TO  THAT  SUPPLIED  TO  AFWL  WITH  FINAL  REPORT  ON 
DISTENOEO  MATERIAL  MODEL  DEVELOPMENT. 


•  <■ 

*  COPPER- 

-COPPER  FOAM  IMPACT 

*  *  A 

1  NTEOT  - 

11 

NJEDIT  * 

6 

NREZON  • 

1 

NSFPRAT 

3 

0 

2  TEDITS* 

5.000E-07 

1 .OOOE- 06 

1 . 500E-06 

2. OOOfc-06 

2 . 500F-06 

3.000E- 

06 

3.500F-06 

A.000E-06 

A . 500E-06 

5.000E-06 

5. 500E-06 

3  JEDITS- 

31  AA 

57  70 

83  106 

A  NTR  ■ 

1 

5  JREZCN* 

6A 

6  NEOTM  » 

10000 

NED  IT  * 

10000 

NPERN  - 

3 

7  STOPS 

JCYCS  - 

*000 

CKS  - 

3-0 

TS  - 

6.000F-06 

8  NMTRLS* 

A 

MATFL  - 

2 

UZERO  - 

7.350E  OA 

POLYURETHANE  FM 

RHOS  = 

1.19 

NEQST  -  1 

NPOR  -  A 

NYO  * 

1 

NCON  *  0 

EQSTC  * 

7.8 16E*1 0 

EOSTD  * 

1.956E*ll 

EQSTE  * 

3. OOOE*  1 0 

EOSTG 

M 

2.00OE*OO 

EQSTH  - 

0. 250E+00 

EQSTS  ■ 

2.21AE*11 

YOS  * 

1.000F*00 

MU  - 

l.OOOE *09 

YADO  * 

0. 

COSO  * 

2. OOOE  *00 

Cl  * 

0.050E  00 

C2  * 

0 

TENSI  1)«  - 

- 1 . OOOE  *10 

TFNSI2I-  - 

- 1 .OOOE* 10 

TENSI 31*  - 

-l.OOOFMO 

RHOP  * 

0.  186 

0.202 

0.333 

0.625 

1.205 

1.250 

COSQ* 

2.0 

2.0 

2.0 

2.0 

2.0 

Cl- 

0.05 

0.05 

0.05 

0.05 

0.05 

1  P2  * 

2.800E*07 

GELP  - 

6. OOOE *06 

YADD  - 

0. 

2  P2  * 

3.350E*07 

DELP  *  - 

-l.A00E*06 

YADD  - 

0. 

3  P2  - 

3.A50E*07 

DELP  * 

2. 500E+05 

YADD  - 

0. 

A  P2  - 

8. 000E*08 

OELP  «  -1.900E*0B 

YADD  » 

0. 

EMELT  * 

3.000E*09 

9.000E*08 

0. 

0.5 

0. 

AK  - 

l.000E*09 

MUP  - 

1.000E+09 

YO  * 

0. 

NZONES*  1, 

10 

CELLS  IN 

6. OOOE-Ol 

CM, 

OELFIN* 

A.  000F-02 

COPPER  FLYER 

RHO  - 

8.9A 

NEQST*  1 

NPOR  *  0 

NYO  * 

1 

NCON*  0 

EQSTC* 

I.A32E  12 

EQSTD  - 

2.A63E  12 

EQSTE  * 

5.310E  10 

FQSTG  - 

2.0A 

EQSTH- 

0.25 

EOSTS  - 

1.593E  12 

YO  * 

1.060E  09 

MU  * 

A.570E  ll 

YADD  - 

0. 

COSQ  * 

A.O 

Cl  - 

0.  1 

C2  - 

0. 

TENSI 11-  - 

l.OOOF+ll 

TENS  121- 

0. 

TENS  1 3 1*  - 

- 1 .OOOE* 10 

NZONES*  1. 

10 

CELLS  IN 

3. 250E-02 

CM 

COPPER  FOAM 

RHOS  - 

8.9A0 

NEQST-  l 

NPOR*  A 

NYO  * 

1 

NCON*  0 

EQSTC* 

1.A32E  12 

EQSTD  - 

2.A63E  12 

EQSTE  - 

5.310E  10 

EOSTG  - 

2.0A 

EQSTH* 

0.25 

EQSTS  - 

1.593E  12 

YOS  - 

1.06E+09 

MU  - 

A.570E  11 

YADD  » 

0. 

COSQ  - 

A.O 

Cl  - 

0.1 

C2  * 

0. 

TENSI II*  - 

•l.OOOE  ll 

TENSI 21-  - 

l.OOOE  10 

TENS (31*  - 

■l.OOOE  10 

RHOP  - 

6.530 

6.539 

8.333 

8.75 

9.063A 

9.2 

COSQ  « 

A.O 

A.O 

A.O 

A.O 

A.O 

Cl  - 

0.  1 

0.  1 

0.1 

0.1 

0.1 

l  P2  * 

l.OOOE  09 

OELP  - 

0. 

YADD  - 

0. 

2  P2  * 

6.000E09  DELP  * 

■  1.000E09  YADD  - 

0. 

3  P2  * 

1.200E  10 

DELP  * 

•5.000C  08 

YADD* 

0. 

A  P2  * 

2. OOOE  10 

OELP  - 

•5. OOOE  08 

YADD  - 

7. OOOE  OR 

EMELT  - 

A.500E  09 

2. OOOE  09 

0. 

0.3 

0. 

AK  - 

7.AA0E*ll 

MUP  - 

3. OOOE  +  11 

YO  - 

0. 

NZONES*  1« 

80 

CELLS  IN 

3. 1A0E-01 

CM, DEL  X- 

2.200E-03 

C7  (MN  GAGES 

RHO  * 

1. 190E*00 

NEQST  -  1 

NPnR  -  0 

NYO  - 

0 

NCON  *  0 

EQSTC  * 

7.8  16E*10 

EQSTD  - 

1 . 956E ♦  10 

EQSTE  - 

3.000E*10 

FQSTG  ’ 

m 

2.000E*00 

EQSTH  * 

2.500E-01 

EQSTS  - 

2.21AE+11 

COSQ  - 

2. OOOE 400 

Cl  * 

0.050E  00 

C2  - 

0 

TENS  ID*  - 

3.000E*09 

TENSI 21- 

0. 

TENSI3S* 

0. 

NZONES*  l. 

30 

CELLS  IN 

7 . 500E-01 

CM,  DELX- 

1 . 200E-02 

•  END* 


NOATE-  10/21/68  IDENT-1041  COPPER-COPPER  FOAM  IMPACT  13404 

COMPUTATIONS  WERE  MADE  WITH  SRI  PUFF  1  ON  THE  COC  3200  AT  SRI.  THF  VFRSWN 
STRESS  HISTORY  AT  INTERFACES  BETWEEN  MATERIALS  1  AND  2  IS  FOP  COMPARISON  WITH 
GAGE  RECORD  OF  OAVE  SCHMIDT  FROM  SHOT  13404 

OF  SRI  PUFF  1  IS  EQUIVALENT  TO  THAT  SUPPLIED  TO  AFWL  WITH  FINAL  RFPORT  ON 
OISTENOEO  MATERIAL  MOOEL  DEVELOPMENT. 


*  * 

•  COPPER- 

-COPPER  FOAM  IMPACT 

*  ♦  * 

1  NTEOT  - 

9 

N JED  IT  * 

6 

NREZON  - 

1 

NSFPRAT 

= 

0 

2  TEDITS- 

9.000E-07 

1.000E-06 

1.900E-06 

2.000E-06 

2.500F-06 

3  JEDITS- 

30  40 

90  60 

70  82 

4  NTR  - 

1 

9  JREZCN- 

64 

6  NEDTM  - 

lOO'JO 

NEDIT  - 

10000 

NPERN  * 

3 

7  STOPS 

JCYCS  * 

1000 

CKS  * 

3.0 

TS  - 

3. 500F-06 

8  NMTRLS- 

4 

MATFL  - 

2 

U2ER0  > 

7.430F  04 

POLYURETHANE  FM 

RHOS  * 

1.19 

NEQST  *  1 

NPOR  -  4 

NYO  = 

1 

NCON  -  0 

EQSTC  - 

7.816EM0 

EQSTD  * 

1.996E*U 

EQSTE  * 

3. OOOE* 1 0 

EOSTG  = 

7. OOOF+OO 

EQSTH  - 

0. 250E+00 

EOSTS  * 

2.214E+11 

YOS  - 

1.0006*00 

MU  - 

1.000E*09 

Y ADO  * 

0. 

COSQ  - 

2.000E+00 

Cl  * 

0.090E  00 

C2  - 

0 

TENStll-  — 1.  OOOE *10 

TENSI2)*  -K000E  +  10 

TENSI3)*  ■ 

-1. OOOE  *1 0 

RHOP  - 

0.189 

0.202 

0.333 

0.625 

1.205 

1.250 

COSQ* 

2.0 

2.0 

2.0 

2.0 

2.0 

Cl- 

0.09 

0.09 

0.09 

0.05 

0.05 

1  P2  - 

2. 800E*07 

OELP  - 

6. 000E*06 

YADD  - 

0. 

2  P2  * 

3.350E*07 

OELP  *  - 

-1 .400E  *06 

YAOD  - 

0. 

3  P2  * 

3.4506*07 

OELP  - 

2. 900E*05 

YADD  * 

0. 

4  P2  * 

8.0006*06 

OELP  -  -1.900E*08 

YADO  - 

0. 

EMELT  « 

3.0006*09 

9. OOOE *08 

0. 

0.5 

0. 

AK  * 

l.OOOE *09 

MUP  " 

l.OOOE *09 

YO  * 

0. 

NZONES-  1, 

10 

CELLS  IN 

6. 000E-01 

CM, 

DFLFIN- 

4.000F-02 

COPPER  FLYER 

RHO  * 

8.94 

NEQST-  l 

NPOR  *  0 

NVO  - 

1 

NCON-  0 

EQSTC* 

1.432E  12 

EQSTD  * 

2.463E  12 

EQSTE  - 

5.310E  10 

EOSTG  * 

2.04 

EQSTH* 

0.29 

EQSTS  » 

1.593E  12 

VO  * 

1.060E  0 

MU  - 

4.570E  11 

YADO  « 

0. 

COSQ  - 

4.0 

Cl  - 

0.1 

C2  - 

0. 

TENS I 11*  "1. OOOE* II 

TENSI2I* 

0. 

TENS<3)-  - 

-1.000E+10 

NZONES"  1« 

10 

CELLS  IN 

3. 760E-02 

CM 

COPPER  FOAM 

RHOS  ■ 

8.940 

NEQST-  l 

NPOR-  4 

NYO  - 

1 

NCON-  0 

EQSTC* 

1.432E  12 

EQSTD  * 

2.463E  12 

EQSTE  - 

5.310E  10 

EOSTG  » 

7.04 

EQSTH- 

0.29 

EQSTS  * 

1.993E  12 

YOS  - 

1.06E*09 

MU  - 

4.S70E  ll 

YADD  « 

0. 

COSQ  * 

4.0 

Cl  * 

0.1 

C2  - 

0. 

TENSI1I-  — l.OOOE  ll 

TENSI2J*  - 

■l.OOOE  10 

TENS! 3 )*  - 

-l.OOOE  10 

RHOP  - 

6.430 

6.449 

8.333 

8.75 

9.0634 

9.2 

COSQ  • 

4.0 

4.0 

4.0 

4.0 

4.0 

Cl  « 

0. 

0. 

0. 

0. 

0. 

l  P2  - 

l.OOOE  09 

OELP  - 

1.667E  08 

YADD  - 

0. 

2  P2  - 

6. OOOE  09 

DELP  *  -l.OOOE  09 

YAOO  - 

0. 

3  P2  * 

1.200E  10 

OELP  * 

■9. OOOE  08 

YADD- 

0. 

4  P2  * 

2. OOOE  10 

OELP  - 

0. 

YADD  ■ 

7. OOOE  08 

EMELT  - 

4.900E  09 

2. OOOE  09 

0.03 

0.65 

0. 

AK  - 

7.440E+11 

MUP  * 

3. OOOE* 1 1 

YO  - 

0. 

NZONES-  1, 

96 

CELLS  IN 

1.910E-01 

CMtDELX- 

2.200E-03 

CT  (MN  GAGE) 

RHO  * 

1 . 190E  *00 

NEQST  -  1 

NPOR  =  0 

NYO  - 

0 

NCON  -  0 

EQSTC  - 

7. 8 16E*10 

EQSTD  - 

1 . 996E* 10 

EQSTE  = 

3. OOOE  *1 0 

EOSTG  = 

2. OOOF+OO 

EQSTH  * 

2.900E-01 

EQSTS  - 

2.214E*11 

COSQ  - 

2. OOOE *00 

Cl  * 

0.090E  00 

C2  * 

0 

TENSI1I-  - 

3.0006*09 

TENS! 2)* 

0. 

TENSI3J- 

0. 

NZONES*  1. 

30 

CELLS  IN 

7.S00E-01 

CM,  DELX- 

1 . 200E-02 

♦  END* 


84 


NDATE*  10/21/68  IDENT-1042  COPPER-COPPER  FOAM  IMPACT  11417 

STRESS  HISTORY  AT  INTERFACES  BETWEEN  MATERIALS  1  AND  2  IS  FOR  COMPARISON  WITH 
GAGF  RECORD  OF  OAVE  SCHMIDT  FROM  SHOT  13417 

COMPUTATIONS  WERE  N.'OE  WITH  SRI  PUFF  1  ON  THE  CDC  3200  AT  SRI.  THE  VFRSION 
OF  SRI  PUFF  1  IS  EQUIVALENT  TO  THAT  SUPPLIED  TO  AFWL  WITH  FINAL  REPORT  ON 
DISTENDED  MATERIAL  MODEL  DEVELOPMENT. 

*  *  *  COPPER-COPPER  FOAM  IMPACT  *  *  • 


1  NTEDT  - 

5 

NJEOIT  ■ 

6 

NREZON  * 

1 

NSEPRAT 

m 

0 

2  TEDITS* 

5.000E-07 

1.000E-06 

1.500E-06 

2.000E-06 

2.500F-06 

:  JED1TS- 

60  70 

80  90 

100  106 

4  NTR  * 

1 

a  JREZON- 

94 

6  NEDTM  * 

10000 

NEOIT  * 

10000 

NPERN  - 

3 

7  STOPS 

JCYCS  - 

1000 

CKS  * 

3.0 

TS  - 

3.000E-06 

8  NNTRLS* 

4 

MATFL  - 

2 

UZERO  * 

7.420E  04 

POLYURETHANE  FM 

RHOS  * 

1.19 

NEQST  -  1 

NPOR  *  4 

NYO  - 

l 

NCON  »  0 

EOSTC  - 

7.8 16E* 10 

EQSTO  - 

1.956E  +  U 

EQSTE  * 

3.000E*10 

EQSTC  - 

2.000F*00 

EOSTH  * 

0.250E*00 

EQSTS  - 

2.214E+11 

YOS  * 

l.OOOE+OO 

MU  * 

1.000E*09 

YADO  - 

0. 

COSQ  » 

2.000E+00 

Cl  * 

0.050E  00 

C2  * 

0 

TENSI 11*  - 

-1.000E+10 

TENSI2J*  ■ 

-l.OOOE* 10 

TENS  1 3 )»  ■ 

- 1 . OOOE+IO 

RHOP  » 

0.186 

0.202 

0.333 

0.625 

1.205 

1.250 

COSQ* 

2.0 

2.0 

2.0 

2.0 

2.0 

Cl* 

0.05 

0.05 

0.05 

0.05 

0.05 

1  P2  - 

2. 800E*07 

DELP  - 

6.000E+06 

YAOD  - 

0. 

2  P2  * 

3. 350E*07 

DELP  -  -1.400E*06 

YAOO  * 

0. 

3  P2  * 

3.450E*07 

DELP  * 

2. 500E+05 

YADO  - 

0. 

4  P2  - 

8.000E+08 

DELP  -  -1.900E*08 

YAOD  * 

0. 

EMELT  - 

3.000E*09 

9.000E+08 

0. 

0.5 

0. 

AK  * 

1.000E*09 

MUP  - 

1.000E+09 

YO  - 

0. 

NZONES*  1, 

10 

CELLS  IN 

6.000E-01 

CM, 

DFLFIN- 

4. 000F-02 

COPPER  FLYER 

RHC  * 

8.94 

NEQST-  1 

NPOR  *  0 

NYO  * 

l 

NCON*  0 

EQSTC* 

1.432E  12 

EQSTO  - 

2.463E  12 

EQSTE  - 

5.310E  10 

EOSTG  - 

2.04 

EQSTH* 

0.25 

EQSTS  - 

1.593E  12 

YO  * 

1.060E  09 

MU  * 

4.570E  ll 

YADO  * 

0. 

COSQ  * 

4.0 

Cl  - 

0.1 

C2  - 

0. 

TENSI 11*  -l.OOOE*!  1 

TENSI 2)* 

0. 

TENS  1 3 )*  - 1 .OOOF* 10 

NZONES*  1, 

40 

CELLS  IN 

1.651E-01 

COPPER  FOAM 

RHOS  * 

8.940 

NEQST-  1 

NPOR*  4 

NYO  * 

l 

NCON-  0 

EQSTC* 

1.432E  12 

EQSTO  * 

2.463E  12 

EQSTE  * 

5.310E  10 

EOSTG  * 

2.04 

EQSTH* 

0.25 

EQSTS  * 

1.593E  12 

YOS  * 

1.06E*09 

MU  * 

4.570E  11 

YADO  * 

0. 

COSQ  * 

4.0 

Cl  * 

0.5 

C2  * 

0. 

TENSI II*  - 

l.OOOE  11 

TENS I 2)*  -l.OOOE  10 

TENS  13  I-  - 

■l.OOOE  10 

RHOP  * 

6.543 

6.558 

8.333 

8.75 

9.0634 

9.2 

COSQ  * 

4.0 

8.0 

8.0 

16.0 

16.0 

Cl  - 

0.5 

0.5 

0.5 

0.5 

0.5 

1  P2  * 

l.OOOE  09 

DELP  - 

1.667E  08 

YADD  - 

0. 

2  P2  * 

6.000E  09 

OELP  - 

■l.OOOE  09 

YADO  * 

0. 

3  P2  * 

1.200E  10 

OFLP  -  -5.000E  08 

YADO* 

0. 

4  P2  * 

2.000E  10 

HELP  * 

0. 

YADO  - 

7.000F  08 

EMELT  * 

4.500E  09 

2.000E  09 

0.08 

0.65 

0. 

AK  * 

7.440E* 1 1 

MUP  * 

3.000E+11 

YO  * 

0. 

NZONES*  1* 

50 

CELLS  IN 

1 . 590E-01 

CM, DEL  X- 

2.000E-03 

C7  (MN  GAGE) 

RHO  * 

1. 190E*00 

NEQST  *  1 

NPOR  *  0 

NYO  * 

0 

NCON  *  0 

EQSTC  * 

7.816E*10 

EQSTO  - 

l.956E*10 

EQSTE  * 

3.000E+10 

EOSTG  * 

2.000F*00 

EQSTH  - 

2.500E-01 

EQSTS  - 

2.214E+11 

COSQ  * 

2.000E+00 

Cl  - 

0.050E  00 

C2  - 

0 

TENSI 1 ) *  - 

3.000E*09 

TENSI2)* 

0. 

TENSI 3  )* 

0. 

NZONES*  l, 

30 

CELLS  IN 

7.500E-01 

CM,  OELX* 

1 . 200E-02 

•END* 


8£ 


NOATE-  t 0/25/66 


IDENT-1047 


TUNGSTEN-TUNGSTEN  FOAM  IMPACT 


1340? 


STRESS  HISTORY  \T  INTERFACES  BETWEEN  MATERIALS  1  AND  2  IS  FOR  COMPARISON  WITH 
GAGE  RECORD  OF  DAVE  SCHMIDT  FROM  SHOT  13402 

CONFUTATIONS  WERE  MAOE  WITH  SRI  PUFF  l  ON  THE  CDC  3200  AT  SRI.  THE  VERSION 
OF  SRI  FUFF  I  IS  EQUIVALENT  TO  THAT  SUPPLIED  TO  AFWL  WITH  FINAL  REPORT  ON 


DISTENDED 

MATERIAL  I 

MODEL  DEVELOPMENT. 

•  •  i 

»  TUNGSTEN-TUNGSTEN 

FOAM  IMPACT  •»  *  * 

1  NT EOT  - 

7 

NJEDIT  * 

6 

NREZON  • 

1 

NSFPRAT  - 

0 

2  TEOITS- 

5.000E-07 

1.000E-06 

1 . 500E-06 

2.000E-06 

2.500E-06 

3.000E-06 

3. 500E-06 

3  JEDITS- 

28  34 

40  45 

50  56 

4  NTR  - 

1 

5  JREZON- 

50 

6  NEDTN  - 

10000 

NEDIT  - 

10000 

NPERN  - 

1 

7  STOPS 

JCYCS- 

2000 

CKS- 

3.0 

TS  - 

3.000E-06 

•  NNTRLS- 

4 

MATFL- 

2 

UZERO  - 

7.3UOE  04 

POLYURETHANE  FM 

RHOS  - 

1.19 

NEQST  *  1 

NPOR  -  4 

Nv0  -  1 

NCON  -  0 

EQSTC  - 

7.8 16E+10 

EOSTD  * 

1.956E+11 

EQSTE  - 

3.000E* 10 

EQSTG  - 

2.000E*00 

EQSTH  - 

0.250E*00 

EQSYS  - 

2.214E+11 

YOS  - 

1 .OOOE+OO 

MU  - 

1.000E+09 

YADD  - 

0. 

COSO  - 

2. OOOE+OO 

Cl  - 

0.050E  00 

C2  - 

0 

TENS! 11* 

-l.Or<1E*10 

TENSI2I-  - 1. OOOE* 10 

TENSI31-  ■ 

•l.OOOEMO 

RHOP  - 

0."  V. 

0.202 

0.333 

0.625 

1.205 

1.250 

COSQ- 

2.0 

2.0 

2.0 

2.0 

2.0 

Cl- 

0.05 

0.05 

0.05 

C.05 

0.05 

l  P2  - 

2.800E+07 

OELP  • 

6.000E+06 

YAOO  « 

0. 

2  P2  - 

3.350E*07 

OELP  -  -1.400E*06 

YADO  - 

0. 

3  P2  - 

3.450E+07 

OELP  » 

2. 500E*05 

YADO  ■ 

0. 

4  P2  - 

8.000E+08 

DELP  «  - 

-1.900E+08 

YADD  - 

0. 

ENELT  - 

3.000E*09 

9.000E+08 

0. 

0.5 

0. 

AK  * 

1.000E+09 

MUP  - 

i.ooor.4-09 

YO  - 

0. 

N20NES-  1 

.  10 

CELLS  IN 

6.000E-01 

CM, 

DEL FIN- 

4.007E-02 

TUNGSTEN 

FLYER 

RHO» 

1.960E  01 

NEQST-  l 

NPOR  -  0 

Ny0  -  1 

NCON*  0 

EOSTC- 

3.101E  12 

EQSTO  * 

3.487E  12 

EQSTE  > 

4.620E  10 

EQSTG  - 

1.6? 

eqsth  * 

0.25 

EQSTS  » 

4.265E  11 

YO  • 

P+ 

• 

• 

P"* 

rvi 

m 

o 

MU  - 

1.550E  12 

YAOD  - 

0. 

COSQ- 

2.0 

Cl  - 

0.1 

C2  - 

0. 

TENSIll- 

-2.000E  ll 

TENS! 2 1* 

0. 

TENSI3J- 

0. 

NZONES-  1 

.  10 

CELLS  IN 

5.050E-02 

TUNGSTEN 

FOAN 

RHOS  - 

19.6 

NEQST-  l 

NPOR-  4 

NYO-  1 

NCON-  0 

EQSTC* 

3.101E  12 

EOSTO  - 

3.487E  12 

EQSTE  - 

4.620E  10 

EQSTG  - 

1.62 

EQSTH  « 

0.25 

EQSTS  * 

4.265E  11 

YOS  - 

1.8  12EM0 

MU  - 

1.550E  12 

YAOD  - 

0. 

COSQ  « 

4.0 

Cl  * 

0.1 

C2  - 

0. 

TENS 111- 

-2.000E  11 

TENS! 21*  - 

-5.000E  10 

TENS! 31-  - 

■2.000E  11 

RHOP  - 

14.56 

14.628 

18.45 

18.90 

20.83 

21.0 

COSQ  « 

4.0 

4.0 

4.0 

4.0 

4.0 

Cl  - 

0.2 

0.2 

0.2 

0.2 

0.2 

1  P2  • 

l.OOOE  10 

OELP  ■ 

0. 

YADO  - 

0. 

2  P2  - 

6.500E  10 

DELP  * 

-8.000E  09 

YADO  - 

0. 

3  PI* 

8.200E  10 

OELP  * 

■2.000E  09 

YADD  - 

0. 

4  P2  • 

2.000E  11 

OELP  * 

0. 

YADO  - 

1.200F  10 

ENELT  « 

6.440E  09 

3.000E  09 

0. 

0.5 

0. 

AK  - 

2.100E+12 

MUP  * 

1 . 050E+ 12 

YO  « 

0. 

NZONES-  1 

,  30 

CELLS  IN 

1.620E-01 

CM,  DELX- 

3 . 000 E— 03 

CT  INN  GAGEI 

RHO  r 

1. 190E+00 

NEQST  -  1 

NPOR  -  0 

NYO  -  0 

NCON  -  0 

EQSTC  - 

7.816E+10 

EQSTO  * 

1.956EM0 

EQSTE  - 

3.000E+10 

EQSTG  - 

2.000F+00 

EQSTH  - 

2.500E-01 

EQSTS  « 

2.214E+11 

COSQ  - 

2.000E *00 

Cl  - 

0.050E  00 

C2  - 

0 

TENSIll- 

-3.000E+09 

TENS! 2 1* 

0. 

TENS! 31- 

0. 

NZONES-  1 

t  30 

CELLS  IN 

7. 500E-01 

CM,  DELX- 

1.200E-02 

♦END* 


ND* TE*  10/23/68  IDENT-1046  TUNGSTEN-TUNGSTEN  FOAM  INFACT  13473 

STRESS  HISTORY  AT  INTERFACES  BETWEEN  MATERIALS  1  ANO  2  IS  FOR  COMPARISON  WITH 
GAGE  RECORD  OF  DAVE  SCHMIDT  FROM  SHOT  13473 

COMPUTATIONS  WERE  MADE  WITH  SRI  PUFF  1  ON  THE  CDC  3200  AT  SRI.  THE  VERSION 
OF  SRI  PUFF  1  IS  EQUIVALENT  TO  THAT  SUPPLIED  TO  AFWL  WITH  FINAL  REPORT  ON 
DISTENDED  MATERIAL  MODEL  DEVELOPMENT. 


•  *  •  TUNGSTEN-TUNGSTEN 


1  NTEDT  * 

7 

N JED IT  * 

6 

2  TEDITS* 

5.000E-07 

1 . 000E-06 

1.500E-06 

3  JEDITS* 

62 

71 

80  80 

08  106 

4  NTR  - 

1 

5  JREZON- 

80 

6  NEOTM  - 

10000 

NEOIT  * 

10000 

7  STOPS 

JCYCS* 

2000 

8  NHTRLS* 

4 

MATFL* 

2 

PLEXIGLAS 

RHO  - 

1.180 

EQSTC  - 

7.816E+10 

EQSTD  - 

1.056E*11 

EQSTH  * 

0.250E+00 

EQSTS  * 

2.214EM1 

COSQ  » 

2. OOOE+OO 

Cl  - 

O.OSOE  00 

TENSI 11*  — l. OOOE+IO 

TENSI 21* 

0. 

NZONES*  1. 

30 

CELLS  IN 

l.OOOE  00 

TUNGSTEN  FLYER 

RHO* 

1.060E  01 

EQSTC* 

3.101E 

12 

EQSTD  - 

3.487E  12 

EQSTH  - 

0.25 

EQSTS  * 

4.265E  11 

YO  * 

1.812E 

10 

MU  * 

1.550E  12 

COSQ* 

2.0 

Cl  - 

0.1 

TENSI 11*  — 2.000E 

11 

TENSI2I* 

0. 

NZONES*  1, 

20 

CELLS  IN 

1.140E— 01 

TUNGSTEN  FOAM 

RHOS  * 

10.6 

EQSTC* 

3. 10  IE 

12 

EQSTD  - 

3.487E  12 

EQSTH  * 

0.25 

EQSTS  - 

4.265E  11 

YOS  - 

1.812E+10 

MU  » 

1.550E  12 

COSQ  * 

4.0 

Cl  * 

0.1 

TENSI 1»*  -2.000E 

11 

TENSI2)*  - 

-5.000E  10 

RHOP  * 

14.26 

14.328 

18.45 

COSQ  * 

4.0 

4.0 

4.0 

Cl  - 

0.2 

0.2 

0.2 

1  P2  * 

l.OOOE 

10 

OELP  - 

0. 

2  P2  * 

6.500E 

10 

DELP  *  -8.000E  00 

3  PI* 

8.200E 

10 

OELP  - 

-2.000E  00 

4  P2  * 

2.000E 

11 

DELP  * 

0. 

EMELT  - 

6.440E 

00 

3.000E  00 

0. 

AK  * 

2. 100E+12 

MUP  - 

1.050E+12 

NZONES*  l. 

50 

CELLS  IN 

3. 100E-01 

C7  IMN  GAGE) 

RHO  - 

1 .  lOOE+OO 

EQSTC  * 

7. 816E+10 

EQSTD  * 

1.056E+10 

EQSTH  * 

2.500E-01 

EQSTS  * 

2.214E+11 

COSQ  * 

2.000E+00 

Cl  - 

0.050E  00 

TENSI 1 )*  - 

3.000E400 

TENSI 2) * 

0. 

NZONES*  1, 

30 

CELLS  IN 

7.500E-01 

FOAM  IMPACT  *  *  * 


NREZON  - 

1 

NSFPRAT 
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